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Machines for Solving Algebraic Equations 


1. Introduction. The search for mechanical means of solving algebraic 
equations has interested mathematicians for well over a century. Two early 
papers date back to the eighteenth century. Perusing a paper of 1758 by 
SEGNER,! in which the author proposes a universal method of discovering real 
roots of equations, based on what we should now call drawing the graph of 
the function y = >-a;x‘, Rown1nG? in 1770 considered the possibility of 
drawing the graph of a polynomial continuously by local motion. Theoreti- 
cally at least, a number of rulers could be linked together so that the pencil 
point on the last ruler would trace the required curve. But mechanical 
limitations of the day caused a reviewer to remark that ‘‘as this is a matter of 
curiosity rather than any use, . . . it is unnecessary to enter any further 
into it at this time.’’ Theoretical methods developed since that day have 
depended for their usefulness on the degree of precision in the mechanisms 
constructed to carry theory into practice, a precision which has greatly in- 
creased in modern times. 

The early mechanical equation-solvers were restricted to finding the 
real roots of equations with real coefficients. But certain electrical methods, 
starting with the one described by Lucas in 1888, were able to handle com- 
plex roots, and even complex coefficients. The modern isograph is an electro- 
mechanical device for finding real or complex roots of algebraic equations. 
In addition to the machines for solving algebraic equations in a single un- 
known, other similar devices have been invented for the solution of simul- 
taneous linear equations in several unknowns. 

Two excellent surveys of earlier mechanisms appeared at the beginning 
of this century, one by MEHMKE® in 1902, revised by p’OcAGNE® in 1909, 
and the other by Moritz‘ in 1905. A few years later GHERs!,' in his book of 
mathematical curiosities, included an illustrated account of some of the 
previously discovered hydrostatic and electric solvers of algebraic equations. 
A comprehensive survey of various dynamical methods of solving algebraic 
equations was given by RIEBESELL® in 1914. The summaries and bibliog- 
raphies published in these papers have been very helpful in the preparation 
of this article, and will be made use of below without further acknowledgment. 

The diverse methods which have been proposed for solving algebraic 
equations mechanically, other than the strictly numerical methods based 
upon the use of calculating machines, fall naturally into about six types, 
and we shall discuss these in the succeeding paragraphs, as follows: (2) 
Graphic and visual methods. (3) Kinematic linkages. (4) Dynamic balances. 
(5) Hydrostatic balances. (6) Electric and electromagnetic methods. (7) 
Methods of harmonic analysis. Of these the first four are usually restricted 
to real roots, whereas the last two may be used to find the complex roots of 
equations. All these types include machines both for algebraic equations in 
one unknown, and for simultaneous linear equations in several unknowns. 

In our description of various devices, it will be less confusing to the reader 
if in most cases we adopt a standard notation for a polynomial whose zeros 
are to be found, which may differ in several instances from those used by the 
authors we quote. Let it be required to determine the roots z = x + ty, 
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(# = — 1), of the algebraic equation 
(1.1) = Co + + + --- + = 0, 


where the coefficients Cn = @m + ib, may be real or complex. The letter R 
will denote any convenient upper bound for the absolute values of the roots 
of f(z). A suitably chosen integral of f(z) will be denoted by F(z), and its 
zeros by Z;, ---, Zn41. To denote real coefficients we shall write a, instead of 
¢m- If only real roots are to be found, the variable z will be called x. Thus 
the notation 


(1.2) f(x) = do + + + --- + a,x" = 0, 


will imply the problem of finding real roots of a polynomial equation with 
real coefficients. In such cases y will often be used to denote f(x). 


2. Graphic and visual methods. Twenty-five years after Rowning’s 
paper, LAGRANGE’ described a graphic method of solving algebraic equations. 
To solve the equation f(x) = 0, Lagrange lays off on the y axis (Lo) the 
n + 1 directed segments OBy = ao, BoB: = a1, Bi Be = az, ---, Bn_1Bn = An. 
The coordinates of the point B,, are seen to be (0, bm), if bm = @o +1 + --- 
+ am. A horizontal line through B, intersects the vertical L;(x = 1) at the 
point C,(1,5,). The line B,1C,, with slope a,, meets a suitably selected 
vertical line LZ, in a point + a,x); a horizontal through 
meets ZL; in Cy_i(1, ba_1 + the line B,2C,_1 with slope a,_1 + a,x 
meets L, in Py_2(x, + + a,x”). Successive points P, are con- 
structed in this way on L, until finally the point P» is found, whose coordi- 
nates are (x, f(x)). The locus of Po, for various lines L,, is the graph of the 
polynomial y = f(x), and the roots are found whenever P, lies on the x-axis. 
Fewer construction lines are involved in the graphic method of LiLL® 
(1867). If we define the algebraic quantities y,, by the successive relations 


(2.1) 0, = — Ym), m=n,n — 1, 2, 1, 
then 
(2.2) = — = + XOn); Jo = Ge — f(x). 


The problem is reduced to constructing successively the segments ym—1, 
m = n,n — 1, ---,1 and finding x by trial so that yo = ao. A rectangular 
framework introduces the coefficients as follows. Starting at O, lay off 
OA, = 4, as a directed segment along the x-axis, lay off AnAn—1 = Gn_1 aS 
a directed segment parallel to the y-axis, lay off An»:An—2 = Gn-2 as a 
directed segment with the positive sense opposite to that of the x-axis; and 
continue at each stage to rotate the positive sense through 90°. Now for 
any assumed value of x draw a line through O with slope — x intersecting 
A,An_. (extended if necessary) at a point P,_: (this is such that An»Pn_1 
will equal y,_1); draw a perpendicular to this line at P,_: intersecting 
Ayn-1An-2 in draw a perpendicular to at etc., until 
finally a point Po is located on A1Ao, such that PoAo = f(x). If x is so chosen 
that Po coincides with A, it will be a root of f(x) = 0. 

This method of Lill was somewhat modified by CREMoNna?® in 1874, and 
was reviewed by Moritz.‘ 
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A graphic method published by CUNYNGHAME” in 1886 gives the real 
roots of an equation of the form x" + Ax + B = 0 to two decimals, if the 
curve y = — x" be first drawn on a suitably large scale. The author’s nota- 
tion is cumbersome, but his idea can be expressed simply by the use of 
coordinates. Let a line of slope A be drawn through (0, B). Then the inter- 
sections of the curve y = — x" with this line y = Ax + B define the real 
roots of the given equation. In the case of the cubic, the author discusses 
the computation of the “impossible” (what we should call imaginary) roots. 
Let P(x, y) be a real intersection of the line y = Ax + B with the cubic 
y = — x’, and let m be the slope of the line drawn from P, tangent to the 
cubic at another point. Then it is easily shown that the three roots of the 
equation are x and —(x/2) + Vm — A. Hence if m < A, two roots are 
imaginary, but these are readily computed from the two slopes. The author 
suggests the use of a special protractor on which slopes can be read directly, 
and marks the abscissas directly on the curve C instead of along the x-axis 
to simplify reading of the roots. 

Extending the graphic idea from the trinomial to numerical equations 
of four or five terms, MEHMKE" describes an apparatus which was displayed 
in a mathematical exhibit in Munich in 1893. The theory is based on the 
fact that if four curved scales in space are cut by a plane, then the four 
readings will satisfy a functional relation. In Mehmke’s model, three of the 
scales are uniform vertical scales A, B, C, parallel to each other but not 
coplanar. A line is determined by a string fastened between a point marked 
u on the A scale and a point marked w on the C scale. A plane is then de- 
termined by viewing this line through a sliding eyepiece set at v on the B 
scale. For an equation of four terms, a non-uniform curved scale, suitably 
graduated, is viewed through the eyepiece and seen to cut this plane in one 
or more points which define the required roots. For an equation of five terms 
a one-parameter family of such curves is used. The curves are constructed 
as follows. Let the given equation be 


(2.3) f(x) = x™ + ux" + vx? + wxt = f. 


On a particular curve corresponding to the parameter f, the point marked 
x is such that its projection on a horizontal plane cutting the A, B, C scales 
in a triangle ABC would be the centroid of masses x*, x”, and x* placed at 
A, B, C, respectively ; and its vertical projection s on one of the three parallel 
scales A or B or C is made to be s = (f — x™)/(x” + x? + x*). Hence the 
intersection of this curve with the visual plane locates the root of the 
equation 
_ + ox? — 
(2.4) x" + x? + x9 


3. Kinematic linkages. The mathematical theory of kinematic linkages, 
discussed long before by Rowning, awaited the day of precision machinery 
before it could be considered practical. In the meantime various theoretical 
devices were discussed in mathematical papers. The equiangular linkage 
described by KEMPE” in 1873 is a device for obtaining real roots of equations 
with real coefficients. It is constructed of m + 1 links Lo, Li, ---, Ln, of 
lengths Jo, 11, ---,7,, joined consecutively at points Po, Pi, ---, P»-1, so that 
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the first link OP» lies along the positive x-axis with its left end O at the 
origin, and so that each of the consecutive links is constrained to make the 
same exterior angle ¢ with its predecessor. The horizontal projection of 
the directed link Ln(= Pm—iPm) is then equal to 1, cos m@, and hence the 


horizontal projection of the directed line OP, joining the two end-points of 
the linkage is 


(3.1) lo + 1, cos @ + 1p cos 26 + --- +1, cos no. 


To find the roots of an algebraic equation f(x) = 0 it is necessary first 
to rewrite the equation in the form (3.1) and then to find those angles ¢; for 
which the free end P, of the linkage lies on the y-axis. If R is any convenient 
upper bound of the roots of f(x), we may set x = R cos ¢, and then calculate 
the Fourier coefficients of the trigonometric polynomial f(R cos ¢). The 
lengths /,, can be taken as any convenient multiples of these coefficients, 
which are linear combinations of the given coefficients a», of f(x). 

Although Kempe did not attempt a description of a mechanical method 
for producing his equiangular linkage, such an attempt was made later by 
BLAKESLEY. In 1907 he published a paper entitled ‘Logarithmic lazytongs 
and lattice-works,”’ and then in 1912 he used the lazytongs idea to devise a 
mechanical construction for Kempe’s equiangular linkage. The construction 
which he gave was erroneous, however, and did not really produce an ap- 
propriate linkage except in special cases. The following is Blakesley’s con- 
struction. On opposite sides of each of the segments P»_1P of Kempe’s open 
polygon OPoP; --- two isosceles triangles PniAmPm and Pm iBmP 
are drawn forming a kite-shaped quadrilateral P»1:AmPmBm. The points Ao 
and By are arbitrary points on the perpendicular bisector of OP», but the 
other vertices A; and B; are chosen successively as the points where A ;_,Pj_; 
and B;_,P;-1, produced, meet the perpendicular bisector of P;_,P;. Then 
2n — 2 rods,—one through A;i1Pi;1:A; and one through B;_,P,_,Bi, 
(i = 1, 2, ---, — 1),—are hinged to each other in pairs at P,;_, and are 
further hinged at their common extremities A; and B;. Two rods connect Ao 
and By, to O, and two rods connect A, and B, to P,, but the vertices of 
Kempe’s polygon are not connected rigidly to each other by rods P;_;P;. 
The joints P; in this configuration can be shown to form an equiangular 
linkage as the lattice is articulated. Blakesley’s error arose in asserting that 
the ratios of the segments OP», PoP:, P:P2, ---, etc., remain constant. If 
this were true we should indeed have a machine for applying Kempe’s 
theory. But in fact, only the ratios of alternate segments remain constant in 
the general case, although all ratios happen to do so in such special cases 
as the logarithmic lazytongs in which the successive lengths form a geometric 
progression. Blakesley’s device fails as an exact solver of algebraic equations. 
It may be of interest nevertheless to note Blakesley’s alternate method of 
transforming the given equation f(x) = 0 into a Fourier series. This he does 
by setting x = tan ¢, and computing the Fourier series for the trigonometric 
polynomial (cos” ¢)-f(tan @). Separate linkages represent the cosine and 
sine terms of the series, and their extremities are fastened at right angles to 
each other to produce a single kinematic linkage. 

An entirely different mechanism for solving equations of the n*» degree, 
based essentially on the graphic method similar to that of Lagrange (1795), 
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was published in the same year (1912) by Mutrueap."* First, by a trans- 
formation of the polynomial which replaces a negative root by a positive 
root, or a root <1 by its reciprocal, the problem is reduced to finding a root 
>1. Two parallel shafts s and s’ are separated by the variable distance 
x — 1. In the same plane, on a line PP’ perpendicular to s and s’ (which cuts s 
in P and s’ in P’, and which we shall call the axis), two centers O and O’ are 
fixed so that OP = P’O’ = 1, OP’ = PO’ = x. To find a root >1 of the 
equation (1.2), mark a directed distance PP, = a, on the shaft s and fix 
it mechanically, let a rod OP, intersect s’ in P;, lay off and fix mechanically 
the directed distance P,P’ ,-1 = @,-1 on the shaft s’, let a rod O’P’,_, 
meet s in P,_;, and lay off Pa_1P»-2 = Gaz on s, etc. To avoid mechanical 
interference, the successive rods can be placed in parallel planes instead of in 
the same plane. The machine is articulated by pulling the shafts s and s’ 
apart. Whenever the point P» or Po lies on the axis, the corresponding root x 
can be read off. 

In the same paper Muirhead described a similar mechanism for solving a 
set of simultaneous linear equations. Each of the variables x, y, z, etc., re- 
quires a pair of shafts (X, X’, and Y, Y’, and Z, Z’ etc.), separated by a unit 
distance, whose distances from a fixed reference shaft are constrained to be 
x and x+1, y and y+1, z and z+1, etc. Each equation 
for example) requires a set of rods of which consecutive pairs are hinged on 
the X, Y, Z shafts, and are forcibly separated by directed distances a;, bj, ¢; 
along the X, Y, Z shafts. The first rod is fixed at O on the reference shaft, 
and the last is separated by the distance d; from O on the reference shaft. 
Ways are devised to avoid mechanical interference. 

Two years earlier (in 1910) NABAUER'® described a rather different 
kinematic machine for solving simultaneous linear equations. The funda- 
mental unit in his machine consists of a horizontal rod bearing two equal 
gears of radius r in vertical planes, one of which rests and rolls on a horizontal 
turntable, and can be slid back and forth on its axis so as to rest at a re- 
quired distance ar from the center of the turntable. As the turntable turns 
through an angle x, the two equal gears turn through an angle ax. To avoid 
slipping, nine gear circles, for a = 1, 2, 3, ---, 9, are provided on the turn- 
table. A separate turntabie mechanism is provided for each of the variables 
x, y, z, and the contributions ax, by, cz are added by differential gears to 
produce a given sum /. Three such sums a,x + by + cz = 1; can be formed 
with three sets of turntables, mounted so that each of the variables x, y, z 
has a vertical axis for itself. If the sums /;, /2, 13 are prescribed, then the 
variables x, y, z are determined. Since friction alone is insufficient to produce 
the desired constraint, gears are used. Since the gears restrict the coefficients 
a;, b;, c; to one-digit integers, a method of successive approximation is out- 
lined which overcomes this difficulty. To a person familiar with early models 
of the differential analyzer, this machine demonstrates some of the same 
ideas and the same difficulties. 

Two decades later a large machine for the mechanical solution of simul- 
taneous equations was built at the Massachusetts Institute of Technology. 
An interesting description of the machine by WILBUR" is accompanied by 
plates showing the machine and some of the details of its operation. Ten 
plates P;, one for each of ten homogeneous variables, are mounted to rotate 
about parallel horizontal axes fastened in a rigid frame. Each plate contains 
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ten slots perpendicular to its axis, nine of which are provided with mid- 
runners on which a positive or negative coefficient a;; may be set approxi- 
mately by hand, and then adjusted with micrometer screws to .0005 in. or 
4 significant figures. The front slot on each plate is used for an accurate 
reading of sin 0@;, the sine of the angle of rotation of the plate from the 
horizontal. The nine slots in a given plate are each associated with one of 
nine given linear equations—which can be considered as homogeneous in 
ten unknowns to be read as sin @;, or as non-homogeneous in nine ratios of 
these unknowns. The ten midrunners corresponding to a single equation lie 
in a vertical plane and are connected by a 60 ft. flexible steel tape. Fastened 
at a point A, the tape rides on 10 ball-bearing top runners, from each of 
which it loops down to a pulley on a midrunner. Then it passes down through 
a clamp C to the bottom, where it passes under a set of 10 bottom runners, 
from each of which it loops up to a second. pulley on the midrunner, and is 
finally made fast at its end B. When the clamp C is released, all the plates 
are free to rotate arbitrarily, since any decrease in length of an upper loop is 
compensated by an increase in the corresponding lower loop, and vice versa. 
The total lengthening of CB, measured along the tape, is 


(3.2) sin 6; 2D; sin 


if the ‘‘coefficients’”’ a;; are set as displacements on nine of the plates P;, 
and the “constants” D; are set as displacements on the tenth plate Po. To 
solve the equations 


(3.3) Lai; D; = © 


the clamp C is made fast and one of the plates is rotated with an oscillating 
motion which constrains the others to move with it. Then the unknowns 
x; = (sin 6;)/sin 0) may be read from the machine. 

Having found x; (usually to within 1% of the largest unknown, but this 
depends upon the stability of the system), the left-hand sides of (3.3) are 
computed accurately to as many figures as desired, with a computing 
machine. If these values d; do not approximate 0 sufficiently, they may be 
reset on the ‘‘constant” plate P» as new constants in a set of equations with 
the same coefficients a;;, of which the solutions are the errors of the first 
approximations. By repeating this procedure, as high a degree of accuracy 
may be obtained as is desired. 


4. Dynamic balances. A paper by BERARD"™ in 1810 was the first to 
suggest using an ordinary balance for finding real roots of algebraic equa- 
tions. His idea was improved by LALANNE” in 1840. The principle of 
the Lalanne balance was to have weights, proportional to the coefficients 
|am|, exert forces placed at directed distances from a fixed reference line 
of +x” if a» is positive, or of —x” if am is negative, so as to produce 
moments @,x" whose sum should be in equilibrium. For convenience in 
staying within a finite portion of the plane, a weight 10|a,| at distance 
+x*/10 can replace a weight |a2| at distance +x’, etc. Balanced on a hori- 
zontal knife-edge is a horizontal rectangle in which guide curves, representing 
+x” divided by a suitable power of 10, are marked either by rigid curved 
wires or by slots in a sheet of metal. From the guide curves hang appro- 
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priate weights representing the coefficients. The constant term ¢o is repre- 
sented by a weight |co| hung on one of two fixed hooks on opposite sides of 
the knife-edge. By means of a horizontal indicator slide, perpendicular to 
the knife-edge, all the hooks are constrained to hang at the same distance 
x from an initial reference line. As the slide is moved from left to right, the 
real roots are read off at positions of equilibrium. 

COLLIGNON,”° in 1873, extended Lalanne’s balance in a theoretical way 
to a complex variable, hanging weights at points representing the complex 
numbers 1, z, 2, etc.; but he was not successful in devising a practical 
machine to move the hooks continuously in the prescribed manner. In 1881 
EXNER* devised a turning balance in which the weights are placed on 
spiral curves. 

Systems of levers were hinged together in a machine designed by Boys” 
in 1886. In this machine a set of m + 1 horizontal axes are placed at levels 
0, 1, 2, --- m, above a fixed pivot. The even-numbered axes are equally 
spaced in a fixed vertical plane s and the odd-numbered ones are equally 
spaced at intermediate levels in a parallel vertical plane s’ at distance 
x + 1 from s. Intersecting the m*" axis at right angles there is a horizontal 
lever with weight pans at distances marked +1 and —1, of which the one 
whose sign is that of a,, is to carry a weight |a,,| so as to produce a moment 
am. Positive moments are clockwise on the left-hand even axes and counter- 
clockwise on the right-hand odd axes. The m—1** lever is supported by a 
sliding joint fastened to the positive pan on the m** lever which is at dis- 
tance x from the m—1** axis. Thus the weight |a,| on the m** lever produces 
a moment a,x on the  — 1*® axis, to which is added the moment @,_; from 
the scale pan. A total moment of (@,_1 + a@,x)x is communicated to the 
n — 2** lever, to which is added the moment a,_2 of the weight |a,_2|, etc. 
Finally, on the bottom 0** lever the total moment is f(x). The two vertical 
planes s and s’ are to be moved apart until equilibrium is reached. Real 
positive roots can thus be read off within limits determined by the size of 
the machine. 

Ten years later (1896) Grant” described a similar machine, with which 
roots from 1 to ~ could be read off on a reciprocal scale, on which the num- 
ber x is marked at a distance 1/x from the mark . Instead of alternating 
the m + 1 horizontal axes in two banks, as Boys did, Grant makes them all 
lie in an inclined plane, displacing them horizontally by worm gears from 
an initial vertical reference plane through the 0-axis in such a way that the 
m-axis is displaced m times as far from reference plane as the 1-axis is dis- 
placed. Scale pans are again hung at distances +1 from the axes, but the 
supports come at distances 1/x from the axes. This machine, like that of 
Boys, is used only to find the real roots, which are read at equilibrium posi- 
tions. It is interesting, in view of later inventions, to note Grant’s statement 
that ‘“‘an imaginary root is as impossible mechanically as it is arithmetically.” 

Although a paper published by Skutscu™ in 1902 refers primarily to 
hydrostatic balances, and will be discussed later in that connection, it 
should be mentioned here for its description of a kinematic circular balance. 
Like spokes of a wheel, a set of m + 1 rods of length r emanate from a center 
O to points Ao, A: --- A»; the whole figure being in a vertical plane. At each 
endpoint A,, is a geared pulley wheel (also in a vertical plane) around which 
is hung a chain suspending a weight pan with a weight adjusted to produce 
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a moment a, on the m* gear. A rod C,,D,, of length 2/ is fastened rigidly to 
this m** gear at its midpoint A,,. Connecting adjacent points C,, and Dn41 is 
a link C,,Dm41, making an angle y with C,,D,, and an angle 6 with DniiCnii. 
The ratio x = sin y/sin 6 is adjusted to be the same for each of the x links. 
Then a balance of moments is obtained when x is a root of f(x) = 0. 
Still another dynamic root-finder was described by PEpprE”* in 1912. 
A set of m + 1 pulleys, numbered from 0 to , are so arranged in a vertical 
plane that the m** pulley is fixed, the 2 — 1** is supported at its center by a 
string passing over the m*® pulley, and each successive pulley is supported 
by a string passing over the previous one, until the string over the last 
pulley (numbered 0) carries a weight to pull the strings all taut. The other 
ends of the m + 1 strings are unwound from drums and are all kept parallel 
at an inclination of @ with the horizontal, where x = 1 + sin @ is the root 
to be found. If a length of string & is unwound from the pulley numbered m, 
then it can be shown that the pulley numbered m — 1 will unwind a length 
of string kx, and that this will finally cause a length kx™ to unwind from the 
bottom pulley (numbered 0). Hence if the machine is set by unwinding 
lengths of string equal to the coefficients @,,@,-1, «++, @ from the drum ends 
of the respective strings, the total displacement of the weight will be f(x). 
At angles @ for which this displacement is 0, the root is read asx = 1 + sin @. 
The method is obviously limited to finding real roots between 0 and 2, but 
can be used for finding the real roots of any polynomial if suitable preliminary 
transformations are first applied to the polynomial. A modification of this 
principle uses springs instead of gravity to pull the pulleys into position, 
and the pulleys are all placed on a movable arm which can turn from @ = 0° 
to @ = 90°. 

For further details on various dynamic balances, the reader is referred 
to Riebesell’s paper.® 


5. Hydrostatic balances. Closely related to the dynamic balances just 
discussed are the hydrostatic balances in which upward forces can be intro- 
duced which are proportional to the volumes displaced by appropriately 
chosen solids. A description by Ghersi of some of these methods has already 
been referred to above.® 

The method proposed by DEMANET” in 1898 is adapted to the solution 
of the cubic equations x* + x = c > 0 and x* — x = ¢ > 0. A cylindrical 
vessel and a conical vessel, both open at the top, are connected underneath 
by a pipe so that water in the two vessels will reach the same level. If the 
tangent of the semi-vertical angle of the cone is ¥3/z, then the volume of 
liquid in the cone will be x* when the surface is x units above the vertex. 
If the cylinder has unit cross-section, it will then contain x cubic units of 
water. Hence if c cubic units of water be added above the level of the vertex 
of the cone, the depth x will be the required root of x? + x = c. For the 
equation x* — x = c, a cylinder is used to displace a volume x of liquid 
within the cone. The cubic y* + py = g can be solved by first writing 
y = 

MESLIN’s”’ contribution in 1900 was a horizontal dynamic balance. 
Solids whose volumes, cut off at a level x, are proportional to x, x’, x*, --:, 
are hung on a lever at directed distances a, @2, a3 from the fulcrum, and the 
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constant d» is represented by a weight |ao| in one of two pans fastened on 
opposite ends of the lever. Before adding the weight |ao|, the weights of the 
various solids have to be balanced by a suitable weight in a scale pan. Then 
the solids are suspended in two interconnected vessels, the weight |ao| is 
added to the proper scale pan, and water is made to flow from a third vessel 
until the balance is again restored to equilibrium. The depth x to which the 
solids are submerged gives the root. In his paper Meslin discusses in some 
detail the exact dimensions of the different solids to be used in his balance. 

Two hydraulic methods were proposed by Emcn** in 1901 for extracting 
the 2** root of any number. The first method, which was extended in a second 
paper to the solution of an arbitrary equation of the m** degree, was es- 
sentially the same as that of Meslin, and involved the depth of immersion 
of certain solids in water. The other, less accurate, depended upon the time 
required to empty a suitable vessel (whose interior is a surface of revolution), 
through a small orifice of area a in the bottom. 

The radius r at height x is obtained as follows: 


(5.1) t= Vx, dt = xrdx/aV2gx = (Vx/nx)dx, 
or 
(5.2) = (aV2g/nx) Vx/Vx. 


A new approach to the same problem was made in 1902 by Sxutscn™ 
who considered not only horizontal positions of equilibrium on a hydrostatic 
balance, but also positions of equilibrium at an angle sin » with the hori- 
zontal. After defining a balance as ‘‘a kinematic chain set up for the applica- 
tion of arbitrary forces, and used for finding relationships between the 
magnitude and position of these forces in equilibrium,’”’ Skutsch describes 
balances of Massavu?® (1878), Grant, and Meslin. He then uses a construc- 
tion, similar to that of Meslin, but instead of the single depth variable x, he 
uses two variables ~, 7, and solves the equation }°a,(t + a,n)’ = 0. By 
means of a float, the inclination sin 7 is related to the depth by the rela- 
tion » = §/c. This makes it possible to get the roots as functions of 7 and 
avoids the difficulties created by negative coefficients. 

A hydrostatic machine for solving a system of simultaneous linear equa- 
tions was developed by Fucus in three papers.***" The last and longest 
of these contains the important elements in the description of his machine. 
A system of stationary containers, partly filled with water, are arranged in 
rows and columns. Each row corresponds to an equation ax+by+---=1;. 
Each column corresponds to a variable x, y, etc. The containers in any 
column are interconnected so that the water will rise in all of them an 
amount equal to the corresponding unknown: x, y, etc. In each row, hollow 
floating cylinders with cross-sections a;, b;, are fastened to a beam which is 
loaded with a weight proportional to /;. The cylinders are thus depressed into 
their containers so that the total displaced water, nearly ax + by + ---, 
is balanced by the load /;. Provision is made for negative coefficients by 
letting water flow into the corresponding float and weight it down, instead 


of being displaced by the float. 


Other variations of the same idea provide for inner and outer water 
levels for each term, one to be filled for positive and the other for negative 
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terms. Changes in level indicate the values of the variables. Another varia- 
tion interchanges the role of rows and columns so that a loading /; of water 
is run into the inner cylinders of one column and the values of the unknowns 
are read off as the weights thus applied to the bars in the several rows. Still 
another machine uses a system of containers in which water communicates 
by rows, and oil, floating on the water, communicates in columns. 

More recently (1940) these hydraulic principles were applied by Scuvu- 
MANN® in designing a mechanism for calculating regression equations and 
multiple correlation coefficients, which is also adaptable to the solution of 
simultaneous linear equations. He estimates that a skilled computer work- 
ing with slide rule and adding machine requires 5m? + m*/4 minutes to 
solve a problem in m variables, whereas his machine will solve the problem 
in 6m minutes. To each independent variable corresponds a horizontal beam 
B; (somewhat like the plates in Wilbur’s machine!’) from which the given 
values y;; of the unknowns y; protrude as directed spikes in a horizontal 
plane. On the end of each spike is a pulley P;;. A rotation of 8; in the beam 
for a variable y; lifts the 7** pulley a distance y;; sin 6; and releases twice 
that length from a cord supporting the pulley. Each equation is represented 
by acord. The 7“ cord, with one end displaced upward by an amount x;, then 
comes down to engage the pulley on the spike yj2, then up over another 
pulley and back down for the next variable ¥;3, etc., and finally ends in a 
heavy float F; suspended in mercury. If we set k; = 2 sin 8;, then the float 
F; will be lifted an amount x; — )cyi;k; when the beams B; rotate through 
8; and the hooks for the dependent variables are raised by x;. The tension 
in each cord is proportional to that displacement. If there are more equations 
than independent variables, the coefficients k; which the machine indicates 


are shown to be the regression coefficients. By articulating one beam at a © 


time with the others clamped, a set of rotations y; can be found from which 
the multiple correlation coefficients can be computed with relative ease. 


6. Electric and electromagnetic devices for finding real and complex 
roots of algebraic equations. A fascinating series of papers published in 1888 
by Lucas gives the reader an opportunity to watch an important new scien- 
tific idea being born. A first paper** presents a generalization of Rolle’s 
theorem. In this paper a set of unit masses, placed at the zeros Z1, Z2, ---,Z» 
of a polynomial F(z) in the complex plane, are assumed to repel a unit mass 
N at z with forces inversely proportional to the distances |Z; — z| from N. 
The components P and Q of the total force are given by the equation 


(6.1) P — iQ = F(z) /F(z). 


Roots of the derived function f(z) = F’(z) are equilibrium positions for N. 
The orthogonal trajectories of the isodynamic lines are curves of degree 
2p — 1 (called “‘lignes halysiques’’) passing alternately through the # roots 
of F(z) and the m = p — 1 roots of F’(z), and are used in Lucas’ generaliza- 
tion of Rolle’s theorem. 

A second paper** considers placing electrodes instead of masses at points 
Z1, Z2, «++, Zn41. The equipotential curves are Cassinian ovals whose nodal 
points occur at the roots of the derivative F’(z), which we denote by f(z). 
In the experimental method due to GuEBHARD,“ a large circular conductor 
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is used in place of a Cassinian oval at the outer boundary of a salt solution 
in which a thin plate of polished metal is immersed. The points of a bundle 
of wires are then placed as electrodes representing as nearly as possible the 
points Z1, Z2, ---, Zn41. Colored rings deposited by electrolytic action fur- 
nish a diagram of equipotential lines whose nodal points are the n roots 
Z1, 22, °**, Zn Of the derivative F’(z). Such figures, obtained experimentally 
for the second degree binomial and the fourth degree trinomial, are displayed 
in the above quoted work by Guébhard.“ 

A third paper of Lucas* gives an electric resolution of algebraic equations 
with numerically given real or complex coefficients, by reducing the given 
equation to an equation of lower degree. An equation f(z) = 0 of degree n, 
can be solved electrically, as described above, if we know the + 1 roots of 
the function F(z), which is an integral of f(z). In applying the first reduction 
method of Lucas we suppose some integral of f(z) to be separated into even 
and odd functions ¢(z*) + 2)(2*), and we then calculate one of the [n/2] 
roots of the polynomial ¥ considered as a function of 2*. If this root be )? 
then a particular integral of f(z), namely F(z) = ¢(z*) — $(A*) + af(2*), is 
divisible by z* — \*, and the quotient is a polynomial of lower degree, m — 1, 
whose roots together with +A and --A are the required electrode points for 
solving f(z). For example, in solving the biquadratic equation f(z) = 5z* 
+ 102° + 322 — 22 — 6 = 0, we have y(2*) = 2* + 2 — 6, which has the 
factor 2? — 2. Then choosing F(z) = (2? — 2)(z* + 32 + 2.52? + 4) we next 
find three roots of the cubic factor. An integral of the cubic is [(5/6)z* + 4] 
X (.322 + 2 + .36), of which the first factor was the “y’”’ function for this 
cubic. Using the four roots of this factored biquadratic as electrode points, 
we solve the cubic electrically ; then we use these three nodal points and the 
two points +V2 as five electrode points to solve the given equation electri- 
cally. An alternative method of reduction is also given, which has the dis- 
advantage of introducing extraneous roots, and which we shall not describe 
here. 

A fourth paper *” discussing the electric determination of the isodynamic 
lines of any polynomial—which are the equimodular lines of its logarithmic 
derivative—describes the field of force due to a system of repulsive centers 
at the roots of F(z), combined with a system of attractive centers at the 
roots of the derivative F’(z). 

The climax of these investigations is reached in a fifth paper,** in which 
the use of partial fractions provides Lucas with a simple and direct pro- 
cedure for locating the complex roots of a polynomial with real coefficients. 
An auxiliary polynomial p(z) = (z — A1)(z — As) — Anys) is chosen 
with distinct roots at convenient points Z on the x-axis. Instead of working 
with f(z) and using the roots of its integral as electrode points, Lucas now 
works with the function f(z)/p(z) = Diu;/(z — A;) and its exponential 
integral 
(6.2) F(z) = — (2 — Angi) 

The constants (residues) w; = f(A;)/P’(A;) are easily computed and are 
introduced into the machine as charges on the fixed electrodes at the points 


z = \;. The roots of the logarithmic derivative of F(z) are the required roots 
of the given polynomial f(z), and appear as nodal points of the isodynamic 
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lines. Experimentally these can be deposited electrolytically on a polished 
metal plate, or they can be traced out by the Kirchhoff and Carey Foster 
method. In the perfected model of the Lucas electric equation-solver, given 
in a sixth paper,*® the need for an outside conductor is eliminated by choos- 
ing n + 2 instead of m + 1 sources and sinks on the axis. The algebraic 
sum of the charges on the m + 2 electrodes is zero, and the apparatus is 
somewhat simpler to construct. Still another paper published two years 
later*® suggests a method of introducing the coefficients into the machine 
by means of resistances proportional to 1/y;. 

An electric device proposed by Kann® in 1902 was the result of his 
inventing a balance similar to that of Lalanne,’* without having been at 
first aware of Lalanne’s result. When Mehmke called this to his attention, 
he conceived the idea of carrying out the same principle electrically, 
replacing the moments of weight by variable resistances. In his mechanical 
balance are slits or wires in the shape of curves +x, +x*/10, +x°/100, etc., 
fastened in a horizontal plane from which are suspended weights proportional 
to the coefficients. These weights are constrained to lie under a movable 
guide slot perpendicular to the x-axis on which the plane is balanced. 
Values of x for which the moments are in equilibrium are the required roots. 
In Kann’s first electrical model, each of the functions +x, +x*/10, etc., is 
represented by a curve of heavy wire fastened in a removable template, 
which is to be mounted vertically by sliding it into one of several pairs of 
slots at the top and bottom of a frame. For positive terms the templates are 
slid in right side up, and for negative terms up side down. A sliding vertical 
frame in a plane perpendicular to all these templates, and perpendicular to 
the horizontal axis in each, serves as a root finder. The coefficients a, are 
introduced by placing in the root finder, where its plane intersects the m*® 
template, a fine wire whose resistance per unit length is proportional to 
|@m| (multiplied if necessary by the power of 10 used in the reduction of the 
template graph). The current then flows through a length of this wire pro- 
portional to the ordinate, and thus through a resistance proportional to the 
given term. These fine wires of considerable resistance make contact with 
the heavy wire curves of the templates on the one hand, and with leads in 
the root finder frame on the other hand, and the successive templates are 
connected in series in such a way that the connecting resistances are either 
negligible or may be assumed constant and be balanced out by the bridge 
resistance which controls the constant term in the equation. As the root 
finder frame is moved in the direction of the x-axes of the templates, the 
roots are found at points where a bridge galvanometer which balances the 
equation reads zero. At these points an adjusted resistance (equal to the 
fixed resistances in the circuit, plus the constant resistances introduced by 
inverting the templates for negative variable terms) just balances a resist- 
ance for the given constant term plus all the fixed and variable resistances 
connected in series in the main circuit. 

To avoid having to use wires of different resistances for each new set of 
coefficients, Kann describes in a final paragraph of the same paper a gear 
mechanism, one to be associated with each of the templates in the frame, 
whereby a resistance proportional to a given term may be unwound from a 
drum in a continuous manner as the root finder is moved along all the given 
power curves simultaneously. For each template there is in the root finder 
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frame a straight vertical shaft with a notch on top which presses against a 
rigid metal curve in the template (such as the curve y = x*/10). The lower 
part of each of these shafts is provided with gear teeth which mesh into an 
exchangeable gear having a number of teeth proportional to the coefficient 
of the corresponding term in the equation. The motion of the shaft, thus 
multiplied by the coefficient |a,,|, is made to turn a drum mounted on a 
horizontal axis, around which wire is wound whose specific resistance takes 
care of the reduction factor of a power of 10 in the template graph. A weight 
hanging on the wire which unwinds from the drum maintains the upward 
pressure on the shaft, and as the root finder moves, a resistance is unwound 
which is jointly proportional to the coefficient (number of teeth) and to the 
power of x (ordinate of shaft top). These resistances are connected in series 
as in the other model, and the roots are found at the zero readings of the 
galvanometer. The method is limited to the determination of real roots, 
however. 

RussEL. & WRIGHT,* in an electric device constructed in 1909, combine 
the principle of slide rule multiplication with addition and subtraction. 
Multiplication is obtained from a thin insulating template in the shape of 
the curve log (y/k) = — x/n, about which a hundred terms of no. 36 in- 
sulated manganin wire is wound, so that the wires are nearly parallel. The 
area in the interval x; = x = x2 is proportional to yi — ye. By adding a 
fixed resistance in series, the total variable resistance is made proportional 
to yi, and the number y; is placed on a logarithmic scale under x; on the 
axis of abscissas. Powers of a variable x can be multiplied into the coefficient 
resistances by setting them on a parallel logarithmic scale. Contact fingers 
using a tangent scale are used to adjust separately the powers of the un- 
known. Then terms are added electrically by combining currents in series 
(or in parallel, using reciprocal logarithmic scales). Finally the real root is 
obtained when combined currents vanish. 

Shortly after the paper of Russell and Wright, RussELL & ALty“ (1909) 
brought magnetism across the electric root-finding trail blazed by Lucas; 
inventing another machine for determining complex roots of equations, and 
stating that the error in its readings would not exceed 1%. In their electro- 
magnetic method ‘‘the horizontal field due to the earth’s magnetism is used 
in an analogous manner to the conducting sheet in Mr. Lucas’ method. A 
drawing board with a slit cut in it, a few pieces of bell-wire, any form of 
‘charm’ compass, ordinary ammeters and rheostats or lamp resistance 
boards, such as are found in every physical laboratory, can be utilized at 
once for the experiment.” 

Using our previous notation, and assuming real coefficients, let 


(6.3) f(z) = + + --- +432 + 
and let = (2 — --+ (2 — Aan), where the real numbers are sub- 
ject to the condition >A; = — @n-1/a,. Then let the partial fraction ex- 


pansion be f(z)/p(2) = @n + — Aj). 

In the horizontal complex plane let the earth’s magnetic flux H be di- 
rected parallel to the imaginary y-axis, so that the real x-axis is perpendicular 
to the magnetic meridian. Let vertical wires through the points (A;, 0) 
carry currents C;, producing a magnetic force C;/5r; at the point s = x + ty, 
where r; = |z — A;|. The components of resultant magnetic force at the 
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point z are 


Ci y Ci x — dj C;/5 
By adjusting the currents C; to equal 5H-y;/a,, respectively, the neutral 
points of zero force will be precisely those for which f(z) = 0. Only n — 1 
ammeters and m — 1 rheostats are required for the apparatus. 

An electric calculating machine for solving simultaneous linear equations 
was described by MALLock* in 1933. Improving upon Mallock’s experi- 
mental model of 1931, this machine for solving ten equations in ten un- 
knowns, to within an error less than 0.1% of the largest root, was constructed 
by the Cambridge Instrument Co. Ltd. The machine will also give a direct 
solution, by least squares, of a set of equations of condition, without forming 
the normal equations. One closed circuit represents each equation, which is 
first transformed algebraically so that all coefficients are less than unity. 
The coefficients are then represented by the relative number of turns of the 
variable coils in a given circuit, each wound about the transformer corre- 
sponding to one of the unknowns. Negative signs are obtained by reversing 
the current in a given coil, and the answers may be read from voltmeters, 
attached to unit coils about the several transformers. Compensators are 
introduced to balance out the effect of resistances in the circuits. Repeated 
approximations may be used to give greater accuracy. 

A recent development in the mechanical solution of algebraic equations 
is the use of rotating and logarithmically expanding parts to represent argu- 
ments and absolute values of complex roots of equations. In an electric 
machine described by Hart & Travis,“ a set of m + 1 coaxially mounted 
generators Go to G, have their rotors rigidly connected together, but the 
stators of m of them are constrained by gears to rotate through @, 26, -- -, 6 
with respect ‘to the one on Gp which is fixed. The alternating voltage 
E cos (wt + k@) on the kt" generator is then multipled by a, in a coefficient 
potentiometer, and by Z* in a modulus potentiometer. The latter is ac- 
tivated by steel tapes wrapped around a spiral cam whose arc is proportional 
to the kt» power of the angle of rotation of the cam shaft. (Only roots of 
modulus <1 are read directly, the others being obtained from the reciprocal 
equation.) The voltages Ea,Z* cos (wt + k@) are then added in series and 
connected to an indicator, the voltage on which may be written: 


(6.5) (E/2) + 
k 

This vanishes, independently of t, if and only if 

(6.6) = 0, or f(Ze*) =0. 


To find the roots Ze the angle 6 is turned fairly rapidly and the modulus 
cam shaft slowly so that Ze* traces a spiral in the unit circle of the complex 
plane. Points of zero reading are the roots, and they can be read within 2% 
in modulus, and 1% in argument. 

The use of a cam to generate a power of a variable appears also in a 
later article by GREEN,*” whose square root extractor is used to change over 
graphically recorded pressure differences into a graph of rates of flow. One 
indicator contacts the edge of a cam whose contour is r = ké?, while the 
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other is driven by a pulley around the same shaft, producing a square root 
mechanically. 


7. Methods of harmonic analysis. The isograph,* ** ®° designed by T. C. 
Fry of the Bell Telephone Laboratories following a suggestion of A. J. 
KEMPNER in 1928, and the S. L. Brown* ® harmonic synthesizer-analyzer, 
are two modern machines for solving algebraic equations. Both solve the 
equation f(re”) = 0 by summing separately the sine and cosine terms in the 
function f(re®) — co, and mapping this function for fixed r as a curve in a 
complex plane. The isograph can handle only real coefficients, whereas the 
Brown analyzer as modified by Brown and Wheeler® can also solve equations 
with complex coefficients. For further details on these machines the reader 
is referred to the reviews which have appeared in MTAC.": ® 


8. Calculating machines. It should not be forgotten that algebraic equa- 
tions of higher degree in one unknown and simultaneous linear equations in 
several unknowns can be solved numerically by standard methods of algebra 
(such as Horner’s or Newton’s method), using a computing machine to save 
labor in the processes of addition, subtraction, multiplication and division. 
An excellent review of the history of Calculating Machines was written by 
BAXANDALL.® One such machine, the Hamann-Automat of the Deutsche 
Telephonwerke und Kabelsindustrie, Aktiengesellschaft, Berlin, described 
by WERKMEISTER,* has the advantage, in solving simultaneous linear equa- 
tion such as = 0, = 0, ---, that a quotient 5; /a; obtained by a 
division can be automatically transferred from the result register to the 
multiplicand register without being copied, and can then be multiplied 
again by each of the successive coefficients a2, a3, --- in the elimination of 
x, between the first two equations. 

A survey article by LILLEy,®* previously reviewed in MTAC, p. 61-62, 
mentions not only calculating machines, but such machines as the BusH 
differential analyzer.** Interesting as this machine is, its primary purpose 
is the solution of differential equations rather than algebraic equations, so 
we shall not describe it here. 

In conclusion the author wishes to express his thanks to R. C. A. for 
his suggestions and help in compiling the literature, and to B. H. Bissinger 
for his assistance in looking up several of the references. In the footnotes 
which follow, those papers which have not been examined are marked in 
the usual manner. 

J. S. Frame 
Michigan State College 
East Lansing, Michigan 


1J. A. DE SEGNER, ‘‘Methodus simplex et universalis, omnes omnium aequationum 
radices detegendi,’’ Akad. Nauk, SSSR. Leningrad, Novi Commentarit Acad. Sc. Imp. 
Petrop., v. 7, 1758, p. 211-226. 

2J. Rowninc, “Directions for making a machine for finding the roots of equations 
universally, with the manner of using it,” R. So. London, Trans., v. 60, 1770, p. 240-256. 

*R. MEHMKE, “‘Numerisches Rechnen,” Encykl. d. Math. Wiss., v. 1, part 2, 1902, 
p. 1067-1073. French ed. by M. p’Ocacne, Encycel. d. Sci. Math., tome 1, v. 4, fasc. 3, 1909, 
p. 339-340, 432-446. 

*R. E. Moritz, “Some physical solutions of the general equation of the m* degree,” 
Annals Math., s. 2, v. 6, 1905, p. 112-126. 
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I. GHERs!, Matematica dilettevole e curiosa, (Manuali Hoepit), Milano 1913, p. 244-249, 
“Meena fisici per la soluzione dei sistemi di equazioni algebriche.” 
*P, RIEBESELL, ‘Uber Z. Math. Phys., v. 63, 1914, p. 256-274. 
7J.L. LAGRANGE, Legons élémentaires sur "les Mathématiques, 1795. English translation: 
Lectures on Elementary Mathematics, Chicago, Open Court, 1898, p. 124-126. 
SE, Lit, ““Résolution graphique des équations numériques d’un degré quelconque a 
une inconnue,” Acad. d. Sci., Paris, Comptes Rendus, v. 65, 1867, p. 854-857; Nouv. Annales 
Math., s. 2, v. 6, i 359-362. 
*L. Cremona, Calcolo Grafico, Turin, 1874. English translation: Graphical Statics, 
Oxford, 1890, p. 70-76. 
10 CUNYNGHAME, “Ona mechanical method of solving quadratic and cubic equations, 
whether the roots be real or impossible,” Phil. Mag., s. 5, v. 21, 1886, p. 260-263. 
11R. MEBMKE, “Uber einen Apparat zur Auflésung Gleichungen mit vier 
oder fiinf Gliedern, ” Z. Math. Phys., v. 43, 1898, p. 338-340. 
ZA.B. KEMPE, “On the solution of equations by mechanical means,”’ Messenger Math., 
2, 1873, p. 51-52. 
1907, BT, ea “Logarithmic lazytongs and lattice-works,”’ Phil. Mag., s. 6, v. 14, 
H. BLAKESLEY, “A kinematic method of finding the roots of a rational integral 
ities AP Mag., s. 6, v. 23, 1912, p. 892-900 
Murrueap, “A mechanism for solving equations of the degree,” Edinburgh 
Math. So., Proc., v. 30, 1912, p. 69-74. 
16M. NABAUER, “Vorrichtung zur Auflésung eines linearen Gleichungssystems,” Z. 
Math. Phys., v. 58, 1910, p. 241-246. 
17J. B. Wicsur, “The mechanical solution of simultaneous equations,” Franklin Inst., 
J., v. 222, 1936, p. 715-724. 
189 J. B. BERaRD, Opuscules Mathématiques, Paris, 1810. 
19L. L. C. Latanne, “Description d’une nouvelle machine a calcul pour résoudre les 
équations numériques des sept premiers degrés,”’ Acad. d. Sci., Paris, Comptes Rendus, v. 11, 
1841, p. 859-860; and report on this machine by A. Cauchy et al., p. 959-961. 
20 F.C. R. COLLIGNON, Traité de Mécanique, v. 2, Statique, Paris, 1873, p. 347-349, 401; 
second ed. 1886; fourth ed. 1903. 
to K. Exner, Uber eine Maschine cur Auflésung hoherer Gleichungen, Progr. Staatsgym- 
nasium, Bezirk, Vienna, 1881. 
2C.V. Boys, “On a machine for solving equations,” Phil. Mag., s. 5, v. 21, 1886, p. 
241-245-+ plate Ill. 
2G. B. Grant, “Machine for solving numerical equations,’’ Amer. Machinist, v. 19, 
1896, Vy 824-826. 
Sxutscu, “Uber Gleichungswagen,” Z. Math. Phys., v. 47, 1902, p. 85-104. 
2% W. Peppte, “A mechanism for the solution of an equation of the ‘ath degree,” Int. 
Cong. Math., Proc., Cambridge, 1912, p. 399-402. 
38 A. DEMANET ,“Résolution hydrostatique de I’équation du troisiéme degré,”” Mathesis, 
v. 18, 1898, p. 81-83. 
7G, MEsLIn, “‘Sur une machine a résoudre les équations,”’ J. de Phys., s. 3, v. 9, 1900, 
p. 339-343. Also Acad. d. Sci. ., Paris, Comptes Rendus, v. 130, 1900, . 888-891. 
28 A. Emcu, “Two hydraulic methods to extract the n* root o any number,” Amer. 
Math. Mo., v. 8, 1901, 10-12, 58-59 
9 J. Massav, “ sur l'intégration graphique et ses applications,’’ Assoc. d. 
Ingénieurs sortis d. Ecoles Spéciales de Gand, Annales, v. 2, 1878, p. 13-15, 203-281 (figs. 
1-65); v. 7, 1884, p. 53-132 (figs. 66-104); v. ‘10, 1886, p. 1-535+12 plates (figs. 105-310). 
Reprinted. “Ap ppendice,” v. 12, 1889, p. 185-45 7+6 plates (figs. 311-350); reprinted, 
Paris, Gauthier Villars, 1890, 264 p. This is a complete list of references to Massau’s work; 
the particular pages discussing the ap tus in question were not identified. The reference 
as by Encycl. d. Sci. Math.', p. 436, to “‘v. 10 (1886/7), p. 58 et planche hors texte fig. 
6”’ is certainly incorrect. 
30K. Fucus, “‘Naherungsweise Elimination durch Mittelwerte,” Archiv. Math. Phys., 
s. 3, v. 17, 1910, p. 103-105; Jahrb. Fort. Math., v. 41, p. 197. 
319 K. Fucus, “Eine Gleichungsmaschine aus kommunizierenden Gefassen,”’ Oesterr. 
Z.f. Vermessungsw., v. 10, 1912, p. 325-329; J. Fort. Math., v. 43, p. 151. 
x Fucus, “Hydrostatische Gleichungsmaschinen,” Z. Math. Phys., v. 63, 1914, p. 
%T, E. W. Scnumann, “The principles of a mechanical method for calculating regres- 
sion equations and multiple correlation coefficients and for the solution of simultaneous 
linear equations,” Phil. Mag., s. 7, v. 29, 1940, p. 258-273. 
*F, Lucas, “Généralisation du théoréme de Rolle,’”’ Acad. d. Sci., Paris, Comptes 
Rendus, v. 106, 1888, p. 121-122. 
SF, Lucas, “Détermination électrique des racines réelles et imaginaires de la dérivée 
d’un polynéme quelconque,” loc. cit., p. 195-197. 
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*F. Lucas, “Résolution yy = des équations age loc. cit., p. 268-270. 
37F, Lucas, “Détermination électrique des lignes isodynamiques polynéme 
quelconque, p. 587-589. 
%F. Lucas, “Résolution immédiate des équations au moyen de I'électricité,”’ loc. cit., 


645-648. 
3° F, Lucas, “Résolution des équations par I’électricité,” loc, cit., 1072-1074. 
oF, Lucas, “Résolution électromagnétique des équations,” . d. Sci., Paris, 
Comptes Rendus, v. 111, 1890, p. 965-967. 
o J. C. Jamin, Cours de Physig do Paris, v. 4, fourth ed. by 
gov, 16 1, p. 170, figs. 2-3. Also P. E UEBHARD, J. d. Physique, s. 2, v. 
af: Kann, “Zur mechanischen Auflésung von crc Eine elektrische Gleichungs- 
Maschine,”’ Z. "Math. Phys., v. 48, 1902, p. 266-27 
Russet, & A. Wriczr, “The. Arthur Wright electrical device for evaluating 
formulae and solving ap age oe * Phil. Mag., s. 6, v. 18, 1909, p. 291-308. 
“A. RussELt & J. N “An electromagnetic par vd y of studying the theory of 
and ‘of any degree,”’ Phil. Mag., s. 6, v. 18, 1909, p. 802-811. 
. R. M. Mattock, “An electrical calculating machine,” R, So. London, Proc., 
v. 1404, 1933, p. 457-483. 

“HC. Hart & I. Travis, ‘‘Mechanical solution of algebraic equations,” Franklin 
Inst., J., v. = 1938, p. 63-72. 
on oH. D. GREEN, el root extractor,” Rev. Sci. Instruments, n.s., v. 11, 1940, p. 

48 R. L. Dretzoxp, “‘The isograph—a mechanical root-finder,”” Bell Laboratories Record, 
v. 16, 1937, p. 130-134. 
4 R. O. MERCNER, “The mechanism of the isograph,”’ loc. cit., p. 135-140. 
8 “Mechanical aids to mathematics: Isograph for the solution of complex polynomials,” 
Electronics, v. 11, Feb., 1938, p. 54. 
L. BROWN, mechanical harmonic synthesizer-analyzer,” Franklin Inst., J., 
v. 228, 1939, p. 675-694; reviewed in MTAC, p. 127. 
L. Brown & L. L. WHEELER, “A mechanical method for graphical solution of 
polynomials,” Franklin Inst., J., v. 231, 1941, p. 223-243; reviewed in MTAC, p. 128. 
53D, BAXANDALL, “Calculating machines,” Encyclopedia Britannica, fourteenth ed., 
v. 4, , 1929, p. 551- 553. Note bibliography at end. 

“Pp. WERKMEISTER, ‘ ‘Die Auflésung eines Systems linearer Gleichungen mit Hilfe der 
Rechenmaschine ‘Hamann-Automat’,” Z. f. Instrumententechnik, v. 51, 1931, p. 490. 
pe... ‘‘Mathematical machines,’”’ Nature, v. 149, 1942, p. 462-465; reviewed in 

“The differential A new machine for solving differential equations,” 
Franklin Inst., J., v. 212, 1931, p. 44 

References to unreviewed items are as follows: 

57G. Rosén, ‘‘Eine elektromechanische ‘Gleichungswage,’” Elektrotechn. Z., v. 50, 
1929, p. 1726-1727. 

58G. Revesst, “Verso soluzioni meccaniche ed elettriche dei sistemi di equazioni 
lineari,” L’ Elettrotecnica, v. 12, 1925, p. 550-553. 


Zeros of Certain Bessel Functions of 
Fractional Order 


The following tables contain the zeros of J,(x) for x < 25, where 
y= + 1/3, +2/3, 41/4, +3/4. These zeros were cbteined by inverse 
interpolation in a thirteen-place manuscript of these functions, computed 
by the NYMTP. The accuracy of the zeros to 10D is guaranteed, and 
the two additional places have a high probability of being correct. 


Ss 

1 2.78088 77239 95 2.00629 96717 90 

2 5.90614 26988 43 5.12306 27427 46 

3 9.04238 36635 83 8.25795 11756 42 

a 12.18134 15289 55 11.39646 76969 87 

5 15.32136 98260 12 14.53629 98843 38 

6 18.46192 72456 89 17.67675 35868 47 

7 21.60278 44489 13 20.81754 94222 32 

8 24.74382 77961 28 (24.740) 23.95855 34952 86 (23.955) 
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s 

1 3.49100 83741 08 1.05850 82594 04 

2 6.65263 55231 22 4.28405 38127 24 

3 9.80161 23591 40 7.44045 44040 05 

4 12.94703 48891 39 10.58817 91486 60 

5 16.09096 95281 99 13.73311 84505 74 

6 19.23414 17604 82 16.87681 75138 75 

7 22.37687 15748 16 (22.384) 20.01985 75839 86 

S 23.16250 59340 75 (23.169) 
j-1/3,0 

1 2.90258 62484 17 1.86635 08588 74 

2 6.03274 70572 66 4.98785 32314 35 

3 9.17050 66694 64 8.12426 53819 40 

4 12.31019 37716 45 11.26351 48254 28 

5 15.45064 89678 17 14.40377 58801 36 

6 18.59148 63361 81 17.54451 06557 21 

7 21.73254 11617 47 20.68550 48061 24 

8 24.87373 14228 06 (24.871) . 23.82665 62470 57 (23.824) 
s 542/80 

1 3.37561 06526 94 1.24304 62596 19 

2 6.53025 59365 13 4.42912 06776 99 

3 9.67658 06352 38 7.57945 84465 30 

4 12.82060 86784 66 10.72474 69244 99 

5 15.96368 38809 06 13.86837 45833 31 

6 19.10627 35045 92 17.01125 45001 33 

7 22.24858 23933 60 (22.253) 20.15373 45371 51 

8 23.29597 58670 60 (23.300) 


For x > 25, ten or more decimal places in the zeros may be obtained 
from the well-known formula (five terms) for the roots of Bessel functions, 
given below; see, for example, G. N. Watson, Treatise on the Theory of 
Bessel Functions, 1922 and 1944, p. 506; the sixth term was supplied by 
W. G. Bick.ey,! Phil. Mag., s. 7, v. 34, 1943, p. 40: 


"3-27-88 
(6949p? — 153855 p?-+1585743u— 6277237) 
105-235. g7 


(u—1)(70197u*— 24 79316y°+480 10494, 
—5120 62548u+20921 63573) 
315-218. 89 


where B = (s + $v — 3)x; uw = 4v*. The first five terms of the above ex- 
pression will yield at least 10 decimals for roots greater than those given here. 
The values of 8 corresponding to the highest roots given here are noted in 
parenthesis; it is apparent that for x close to 25, 8 approximates the root to 
two decimals at least. 

The author desires to express his appreciation of assistance rendered by 
several members of the NYMTP in checking these values. 


MILTON ABRAMOWITZ 


NYMTP 


1 It may be noted that in Bickley’s article the fifth term of the formula has an erroneous 
number 6277327, for 6277237. 
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RECENT MATHEMATICAL TABLES 


169[B].—BAASMTC, Mathematical Tables, Volume IX. Table of Powers 
giving integral powers of integers. Initiated by J. W. L. GLAISHER, extended 
by W. G. Bickiey, C. E. Gwytuer, J. C. P. E. J. TerNoutn. 
Edited by Miller. Cambridge University Press, 1940, xii, 132 p. 21.5 x 
27.9 cm. 15 shillings. 

NYMTP, Table of the First Ten Powers of the Integers from 1 to 1000. New 
York, December 1, 1938 (title page), 1939 (outside cover), 80 p. 21 XK 35 
cm. Reproduced in mimeographed form. Originally listed at $.50; out of 
print. Not to be reprinted. 


These two tables, published within a year of one another, are similar in their purpose, 
to provide computers with an adequate table of the higher powers of the integers. The first, 
however, covers a much greater range than the second, although for powers in which the 
two overlap the second is somewhat easier to use because of the various printing devices 
which were used to condense the longer work. 

The BAASMTC tables originated with a work prepared by J. W. L. GLaIsHER more 
than 70 years ago to which the following reference was made in the 1873 Report of the 
Committee on Mathematical Tables: ‘“‘Mr. J. W. L. Glaisher has had formed in duplicate a 
table giving the first twelve powers of the first thousand numbers, which, after the calcula- 
tion has been made independently a third time, will be stereotyped and published, probably 
in the course of 1873; it is hoped that it will help to make the tabulation of mathematical 
functions less laborious and difficult.’’ But the tables never appeared. 

Intrigued by this reference and also by a reference to these tables made in 1905 by A. J. 
C. CunnincuaM, L. J. Comrie initiated a diligent search for a copy. His efforts were finally 
rewarded by success and in 1935 H. J. WoopDALL presented to the Committee the copy to 
which Cunningham had made reference. This is a proof copy and is probably the only one 
in existence. The introduction to the work under review states that “without this copy it is 
doubtful if the production of the present volume would have been undertaken.” 

With the original Glaisher computation as a nucleus, the scope of the table was ex- 
tended to give x" for the following ranges: 


m= 2(1)12 x = 1(1)1099 
nm = 13(1)20 x = 1(1)299 
nm = 21(1)50 x = 1(1)120. 


Because of the limitations of space, the powers corresponding to m = 28 and 29 were 
curtailed, and for values of m greater than 30, except for m = 40 and m = 50, the powers were 
also curtailed. In these abbreviated values the last figure printed is always the true value 
and is not raised. For convenience in use the number of figures retained is a multiple of five, 
and the counting of the digits is assisted by small figures printed at the top and bottom of 
the columns or in the margins. In the incomplete values never fewer than 21 figures are 
retained; below 100 not more than 25 are given, but in the range to 120 the number may be 
as great as 28. 

“The tables are arranged so that, whenever possible, an opening contains powers of a 
hundred consecutive integers, and so that consecutive powers of a given integer are normally 
on the same or adjacent openings. . . . The complete range of powers for each hundred 
consecutive integers forms a single ‘chapter’ which is completed before the next is started.”’ 

For compactness in printing various devices have been adopted. Thus, we find that 
“powers have in some cases been combined in pairs so as to give an approximately constant 
total number of digits. For example, in the first ‘chapter,’ x® and x*° appear together on 
pages 4 and 5, followed by x!° and x** on pages 6 and 7, and so on, up to pages 24 and 25, 
which contain x!* and x®°. Thus for consecutive powers of, say, 57, we turn forward until we 
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come to x**, then backward (as indicated by a footnote) until x* is reached; another foot- 
note directs us to page 27 for x*!, and from there we proceed forward to x°®.”” 

Unusual care was taken to insure accuracy in the table. For this purpose the following 
formula was used, 


1 


where we abbreviate: 
n(n — 1)(m — 2) --- (n—1) 
(2r + 2)! 


in which B,,; is the (r + 1)th Bernoulli number. The powers were summed in blocks of 50 
and these totals compared with those computed directly from the formula. The figures were 
also read against those given in other existing tables where the ranges overlapped. This 
resulted in the detection of a number of errors in other tables, which have been listed in the 
“Bibliography.’”? No errors were discovered in the calculations of Glaisher, nor in the 
NYMTP tables, which were loaned by A. N. Lowan for this check. 

In addition to the table of powers there is given a table of the binomial coefficients? (*) 
form = 2(1)50,7 = 2(1)10and a table of the numerators and denominators of the coefficients 
@2r41 for values of m = 50. 

The NYMTP table was the first publication of that Project. It was computed from the 
first four powers of the integers given in P. Barlow’s Table of Squares, Cubes. . . . Because 
of unusual precautions taken to avoid errors in the final stencils, it was believed that the 
table was without error. The check by the BAASMTC, which we have mentioned above, 
confirmed this belief. 


= 


H. T. D. 
1 Errors in the following works are listed: 
P. Bartow, New Mathematical Tables . . ., London, 1814 
K. ——- Sieben- und mehrstellige Tafeln der Kreis- und Hyperbelfunktionen . . ., Berlin, 
1 


K. Hayasat, Fiinfstellige Funktionentafeln . . ., Berlin, 1930 

K. Havasu, Tafeln fiir die Differenzenrechnung . . ., Berlin, 1933 

K. Pearson, Tables for Statisticians and Biometricians, Part II, London, 1931 

H. W. WEIGEL, x* + y" = 2"? Die elementare Lésung des Fermat-Problems . . ., Leipzig, 
33 


19 
L. —. Vollsténdige Tafeln der Quadrate aller Zahlen bis 100 009, third ed., Berlin, 
1938. 


2 A complete table of the binomial coefficients (7) for m = 2(1)50, r = 2(1)50, is con- 
tained in J. W. L. GLatlsHER, Messenger Math., s. 2, v. 47, 1917, p. 97-107. 


170[B].— Square Root Divisors, 1941, 2 p., Square Root Multipliers, 1940, 2 p., 
Cube Root Divisors, 1944, 2 p. Publications MM 56, 57, 68 of the Marchant 
Calculating Machine Co., Oakland, California. 21.5 K 28 cm. 


These tables are one-page substitutes for more extensive tables of square and cube 
roots and are supposed to be used in connection with any standard computing machine. 

The table of Square Root Divisors, T. 56, gives 8-place values of 2A# and 2(10A)}# for 
166 integral values of A (at irregular intervals) between 100 and 1000. The table of Square 
Root Multipliers, T. 57, gives for the same values of A the reciprocals of the values in T. 56. 

To use T. 56 to find an approximation to the square root of a given number N, we enter 
the table at the value of A which is nearest to N and take out the value of 2A# (or perhaps 
2(10A)#). Then the approximation 


2A 


is obtained by machine division of N + A by the tabulated entry. Thus the method is based 
on the fact that the arithmetic and geometric means of two nearly equal numbers are very 
nearly equal. 
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T. 57 merely replaces the operation of division by multiplication. Probably the average 
computer will find T. 56 a bit more efficient than T. 57 when his machine has automatic 
division and the quotient can be copied down while the machine is operating. In some cases 
when the square root is in the denominator T. 57 might prove more effective. 

It may be argued that no table for square root is necessary if one has a computing 
machine. Most computers are familiar with a square root process based on the identity 


which is well suited to some machines and which will yield the square root of an occasional 
number sooner than the operator can locate a table. Tables 56 and 57 will be found more 
efficient when a great many scattered square roots to 5S are needed quickly. 

The Table of Cube Root Divisors, T. 68, gives 6-place values of 3A*/*, 3(10A)*/* and 


" 3(100A)"* for 162 values of A between 100 and 1000. It is used in a similar way to obtain 


cube roots by the approximation 


The irregular intervals of the argument A of the three tables are so chosen that one 
application of the process will give the required root correct to at least 5S, no matter what 
value of N is given. On the reverse sides of T. 56 and 68 are examples of the re-application 
of the appropriate formulae a second time to obtain approximations correct to 9S and 10S. 

The method involved in these tables was apparently suggested by comments of L.J.C. 
in his “Introduction” to the third edition of Barlow’s Tables (1930). The general method, 
which is usually superior to interpolation, may be used with any table of square or cube 
roots, or indeed any table of squares or cubes. T. 68, however, is decidedly more convenient 
for cube roots. 

The reviewer has two photostatic copies of a 3-page manuscript table of the same sort 
for 5th roots prepared and supplied by the same company. No doubt this useful table will 
be published in due course, to replace Basic Publication MM 88, of Aug., 1940, 2 p. 

Quite another method for cube and 5th roots has been given by Dederick.! A one-page 
table enables the computer to obtain either root to 10S. 

1L. S. Deperick, “A modified method for cube roots and fifth roots,” Amer. Math. Mo., 


v. 33, 1926, p. 469-472. See also D. H. LEnmer, “‘On the use of the calculating machine for 
cube and fifth roots,” Amer Math. Mo., v. 32, 1925, p. 377-379. In the “table of cube and 


fifth roots” for this latter article are the following slips: ¥i0, for 93193, read 93192; 41000 
for 71705, read 71706; ¥1.02 for 27709, read 22710; V5, for 29662, read 29661. 


171[D].—A. “Mehrkreisige Siebenschaltungen mit ausge- 
glichener Resonanzkurve,” Elek. Nach. Technik, v. 20, Oct. 1943, p. 
247-248. 21.7 X 28 cm. 

There are tables of exact values, and to 6D, of cos moe, and sin oe, 

n=1, 2, 4, »=1, 2; 2 =5, 6, » = 1(1)3; = 8, » = 1(1)4; = 10, 

» = 1(1)5; m = 12, » = 1(1)6. Also exact values, and to 5D, of cos a a 


24 a” 

» = 1(1)6. 

172[E].—NATIONAL DEFENSE RESEARCH CoMMITTEE, Division 6, Tables of 
the Bipolar Transformation, compiled by M.I.T. Underwater Sound 
Laboratory, Report, November 1944, 13 leaves + 3 plates (2 folding). 
21.5 X 28 cm. Printed by the photo-offset process from manuscript. 


50 
ere 
‘his 
the 

N+2A 
the N's = 
(7) 
nts - 
the 
use 
the 
ve, 
q 
lin, 
rig, 
lin, 
on- 
int 
ibe 
for 
are 
56. 
ter 
aps 


4 


358 RECENT MATHEMATICAL TABLES 


These tables are available only to certain Government agencies and 
activities. 
In all transmission line problems the transformation 


tanh [x(a + i8)] = R + iX = pe'# 
is of importance. 
The tables give values of R(SD), X(5D), p(4D), and ¢ (to the nearest 1’) as functions 
of a = 0(.01).2(.02).4(.05).6(.1)1, and 6 = 0(.02).5. 
The following formulae are useful: 


we sinh (27a) x= sin (278) 
cosh (2ra) + cos (2x8) ’ cosh (27a) + cos (2x8) ’ 
R? + [X + cot (248)? = csc* (2x8), 
[R — coth (2%a)}? + X? = csch? (27a). 


The first plate (27.8 X 34.3 cm.) is a graph of tanh [x(a + i8)] = pe**, ¢ = 0(0°.5)90°, 
p = .1(.01)1(.1)10, a = .01(.01).06(.02).2(.05).3(.1).5, 6 = .01(.01).06(.02).2(.05).3(.02)- 
-44(.01).5. 

The second plate (27 X 27.8 cm.) shows contours of the transformation R + iX 
= tanh [x(a + 78)], R = 0(.05)5, X = 0(.05)5, a = .005(.005).1(.01).25, .3(.02).4(.1).6, 
8 = 0(.02).1, .15(.01).4(.05).5. There are three errors in marking the letters a and 8, and 
7 places where I is substituted for X. 

The third plate (27.8 X 19.1 cm.) displays contours for the equation pe‘* 
= sinh [x(a + 78)], p = .1(.005).2(.01)1(.05)2(.1)8, ¢ = 0(2°)90°, a = .02, .05(.05).3(.1).9, 
B = .02, .05(.05).5. 


R 


173[E].—A. Pucuer, ‘“‘Rechteckplatten mit zwei eingespannten Randern,”’ 
Ingenieur-Archiv, v. 14, 1943, p. 250-251. 19.3 K 26.8 cm. 


There are 8 tables, to 5S, for x = .25(.05)3.6, of the following functions: 


sinh rx + ax cosh xx cosh rx + xx sinh rx wx sinh rx 
ax cosh sinh xx cosh rx 


}sinh 2ex — xx’ }sinh 2ex tex’ }sinh + 


174([E].—P. B. Wricut, ‘Resistive attenuator, pad and network, theory and 
design’”’ part 2 of a 4-part paper, Communications, v. 24, Oct. 1944, p. 64, 
66, 68. 19.6 X 27.1 cm. 


Tables for 20 log e*(e = k? > 1) = 0(.01).2(.05).4(.1)4(.5)30(1)60(5)140, 150 of (a) 
sinh @ = (4% — 1)/2k; (b) cosh @ = (k? + 1)/2k; (c) tanh 0; (d) tanh 40 = (k — 1)/(k + 1); 
(e) cosh? 6; (f) sinh* 40 = (k — 1)*/4k; (g) csch @; (h) sech 6; (i) coth @; (j) coth 40; (k) 
sech? @; (1) csch? 46. These tables are to 5-9S. 


175(I].—A. N. Lowan & H. E. Sauzer. “Coefficients for interpolation 
within a square grid in the complex plane.” J. Math. Phys., M.I.T., v. 23, 
1944, p. 156-166. 17.5 & 25.5 cm. 


The usual theory of polynomial interpolation holds in the complex domain as well as 
on the real axis. Thus, given an analytic function f(z) of the complex variable z = x + iy 
and » + 1 arbitrary points zo, 21, ---,2n within the domain of regularity of f(z), there 
will be one, and only one, polynomial ¢n(z) of mth degree which will assume the prescribed 
values f(z.) for k = 0, ---, m. The degree of approximation of ¢,(z) to f(z) can be judged 
from Taylor’s formula. The polynomial ¢,(z) can be represented in many forms, but we are 
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here concerned only with the standard Lagrange form 


P,(z) 
on(z) = 2 P, (zs) f 


where P,(z) = (2 — Zo) --+ (2 — — ( — 

Originally, this formula has been regarded as a purely theoretical tool. However, in the 
real domain tables of the coefficients P,(z)/P,(z:) are now available for the typical case of 
equally spaced z;. With such tables! one has only to perform the m + 1 multiplications by 
f(z), and the products are accumulated in any modern computing machine without clearing. 
Thus, using tables and a modern computing machine, the Lagrange interpolation formula 
will lead to the required value ¢,(z) without any writing and auxiliary computations. In 
the complex plane the situation is different only in that a multiplication of two complex 
numbers requires four ordinary multiplications. 

The present paper provides tables of the quadratic and cubic Lagrange interpolation 
polynomials if the three or four given points are vertices of a typical square: zo, 21 = zo + h, 
21 = 20 + th, 22 = 29 + h + th (h real), and it is desired to interpolate for a point z = zo 
+ Az within the square. In the usual way one puts Az = h (p + ig) so that O=p=1, 
0 =g = 1. We have then to interpolate at the point P = p + ig of the unit square. The 
quadratic and cubic formulae are written in the form 

= + LPP )f(eo) + LE (er) 
and 
oa(z) = + LP(P)f(eo) + Fes) + LP (PF). 


The seven coefficients L®(P) etc. are tabulated (p. 160-166), each in the form of a 
double-entry table with the arguments ~ and g varying from 0 to 1 in steps of .1. The values 
are exact, which means that the coefficients LY: L®(P), L®(P) are given to 3D and 
LP) to 2D, while the coefficients L’(P), k = — 1(1)2, are given to 4D. 

W. FELLER 
Brown University 


1NYMTP, Tables of Lagrangian Interpolation Coefficients, New York, 1944. 


176{I].—H. E. Sauzer, ‘Table of coefficients for inverse interpolation with 
central differences,” J. Math. Phys., M.I1.T., v. 22, 1943, p. 210-224. 
17.5 X 25.5 cm. 


The author has done a genuine service to practical interpolation by tabulating not only 
the coefficients for inverse interpolation by central differences, as given in the article under 
review, but also for a similar tabulation of the coefficients of the formula with advancing 
differences; see the review by D. H. L., MTAC, p. 315 f. The difficult problem of finding 
inverse values from tables which are computed to arguments that are not sufficiently close 
to assure linear interpolation, is thus made much simpler by means of these tables. 

The tables computed by Salzer are based upon the inverse of what is commonly called 
the Laplace-Everett formula for interpolation. If f(x) is the function tabulated for x = xo, 
xo + h, xo + 2h, etc., and if p lies between 0 and 1, then the Laplace-Everett formula 
may be written: 


(1) + ph) = f(x) + paf(x) — + Ex(p) + 
+ + + h)] — --- 
where g = 1 — , &f(x), 5*f(x) etc. are central differences, and 


Ex(p) = — p*)/3!,  Ex(p) = — #*)(2? — p*)/Stetec. 


If we make the further abbreviations: 


m = Ph) — 


5"f(x) 
Af (x) 


Af(x)’ 


1 


J 
| 
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then equation (1) can be written in the form: 
pb = m+ E,(g)d? + Es(p)d:? — Ex(q)dot — Ex(p)dit + ---. 


This formula can now be inverted by means of Lagrange’s formula for the inversion of 
functions of the form p= m+ A¢(p), subject to the condition |A¢(p)| <r, where r= |p—m]|. 
This inversion leads to the following expansion: 


m+ + — m*)d? + A(m)do' + B(m)di* + C(m) (de)? + ---, 
where »C; is the third binomial coefficient and the functions A(m), B(m), C(m), etc. are 
polynomials in m. If one defines a product of the form (d”)*(d%)* as being a term of order 
br + qs, then the inversion is carried out through terms of sixth order. The original calcula- 
tion of the coefficients by H. T. D.! were checked and one printing error detected, the first 
number in the coefficient of the term do*(d,*)? being 1 instead of 2. 

Exclusive of the coefficients of the terms of second order there are 15 coefficients to be 
tabulated, which the author has designated by A(m), B(m), ---, O(m). The first two coeffi- 
cients mC; and m(1 — m*)/6 are to be found in the recently published Tables of Lagrangian 
Interpolation Coefficients of NYMTP (MTAC, p. 314f), and are accordingly omitted by the 
author. The remaining 15 coefficients are computed to 10D. Those of fourth order, namely 
A(m) through E(m) are given over the range m = 0(.001)1 and the author makes use of the 
relationships: 

A(m) = B(i — m), —C(m) = D(i — m), E(m) = — E(i — m), 
to reduce the space occupied by the tables. The coefficients of sixth order are given over the 
range 0(.1)1. 


1H. T. Davis, Tables of the Higher Mathematical Functions, v. 1, Bloomington, Ind., 
1933, p. 82-83. 


177[J].—H. StenzEt, Schallfeld eines Strahlers in einer Mediumschicht 
mit schallweicher und schallharter Begrenzung,” Ann. d. Phys., s. 5, 
v. 43, June, 1943, p. 30,18. 13.8 X 21.5 cm. 
1 2 3 


1 
Table of He) = 


+--- 


for x = [.005(.005).5; 4D]. 

On p. 18 iH{(x) is given, to 4D, for x = .6084, .4156, .4712i, .9024i, 1.266i, 1.610%, 
1.9432, 2.2727, 2.5952, 2.9181, 3.2387, 3.5581, 3.8767, 4.1947, 4.5127, 4.8297, 5.1457, 5.4622, 
5.7942, 6.094:. 


178[L].—P. K. Bose, “On confluent hypergeometric series,” Sankhyd. 
The Indian J. Statistics, v. 6, Feb., 1944, p. 407-412. 22.5 X 29.2 cm. 


For integral and half integral values of a and c the confluent hypergeometric series 
1F;(a; ¢; x) occurs in F-Statistic under non-null hypothesis, in Studentized D?-Statistic and 
in the expression for the multiple correlation coefficient for a particular type of parent 
population. Tables with 5S are given here for a = 2(1)15 and with 4S for a = 16(1)25. 
The values chosen for c are 1, 2, 3, 4 and the values of x are 2, 3, 3.8, 5.7, 7, 8.4, 9, 10.5, 
13.5. A short table is given also to show the agreement between the value calculated from 
the series and the value calculated by recurrence formulae starting from the known values 
fora = 2,3;¢ = 1, 2. 

These new tables include some values of x not considered in previous tables. In those 
of Wess and Arrey! for c = 1(1)7, a = — 3(.5)4 the range with 4S was x = 1(2)6(1)10. 
In Airey’s first table* for 2c = —3(2)3, 2a = —8(1)8 the range with 6S was x =0(.1)1(.2)3(.5)8. 
In the tables of GRAN OLsson® ¢ has the values 1, 2, 3 but a is not an integer, 4D being given. 
In the tables of CHAPPELL‘ c = 1,a = § — k, k = 1(1)10, kx = .1(.1)1.5, 2(1)10 with 4D. 
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The functions for which c is an integer occur in many physical problems some of which 
are mentioned by Webb and Airey.! The Bessel function I,,(z) in which m is an integer or 
half an odd integer is related to these functions, so also are the polynomials of Laguerre and 
Sonin. The case c = 3 occurs in the work of SEN* on stresses in some rotating circular disks 
of varying thickness while the function 1F;(a; » + 2; —2ky) occurs in the work of Lams* 
on atmospheric oscillations, and some calculations were impracticable on account of the lack 
of tables. The polynomials of Hermite correspond to the case in which 2c is 1 or 3. Some 
tables of these functions are already available. The author indicates that further numerical 
results are forthcoming. 

Some numerical results relating to the roots of Laguerre polynomials L.(x) = F(—n; 1; x) 
suitable for numerical integration over the range (0, ~) have been given by KosHLiaKov.’ 
Tables of the function 


H(m, a, x) = exp (—ix)F(m + 1 — ia; 2m + 2; 2ix), 


with 6D, have been published by Lowan & Horenstetn® for m = 0(1)3, a = 0(1)10, and 
x = 0(1)10. The 2 in 2ix has been omitted on the cover of the reprints. 
H. B. 


1H. A. Wess and J. R. Arrey, “The practical importance of the confluent hypergeomet- 
ric function,” Phil. Mag., s. 6, v. 36, 1918, p. 129-141. 

2]. R. Arrey, “The confluent hypergeometric function,” B.A.A.S., Report, 1926, p. 
276-294; 1927, p. 220-244. 

?R.’ GRAN Otsson, (a) “Biegung kreisférmiger Platten von radial veranderlicher 
Dicke,” (b) ‘“‘Tabellen der konfluenten hypergeometrischen Funktion erster und zweiter 
Art,” Py yore v. 8, 1937, p. 81-98, 99-103. 

‘ CHAPPELL, “The properties of a new orthogonal function associated with the 
confluent hypergeometric function,” Edinburgh Math. So., Proc., v. 43, 1925, p. 117-130. 

5B. Sen, “Note on the stresses in some rotating circular disks of varying thickness,” 
Phil. Mag., s. 7, v. 19, 1935, p. 1121-1125. 

*H. AMB, “On atmospheric oscillations,” R. So. London, Proc., v. 84A, 1910, p. 551-572. 

tN. S. KosHLIAKov, “O vychislenii po formule mekhanicheskikh kvadratur opre- 
delennykh integralov s beskonechnymi pol. w Thee [On the calculation of in s to in- 
finite limits by means of formulae of mechanical quadratures], Akad. Nauk, i ‘ 
Izvestita, s. 7, Fiziko-matematicheskoe otdelenie, v v. 7, 1933, p. 801-808. 

*A.N. Lowan &W. HorRENSTEIN, “‘On the function H m, a, x) =exp (—ix) F(m+1—i0; 
2m+2; 2ix),”” J. Math. Phys., M.1.T., v. 21, 1942, p. 273-283. 


179[L]—BAASMTC, Mathematical Tables, Volume VI. Bessel Functions, 
Part I, Functions of Orders Zero and Unity. Cambridge, University Press, 
1937, xx, 288 p. 21.5 X 28 cm. Compare MTAC, p. 282f. £2. 


The circular cylinder is the most symmetrical body after the sphere. This is perhaps the 
reason why Bessel functions which arise in the solution of boundary problems for domains 
with circular symmetry are the most widely used transcendental functions after the circu- 
lar functions which enter in the solution of boundary problems for domains with spherical 
symmetry. From the standpoint of the frequency of occurrence, the Bessel functions of 
orders zeroand unity are undoubtedly the most important. Moreover, the recurrence formula 
between Bessel functions of three consecutive integral orders make it a relatively simple 
matter to generate the values of Bessel functions of integral orders from those of orders zero 
and unity. These considerations serve to stress the great contribution to science of the tabular 
volume before us, devoted to Bessel functions of orders zero and unity. This volume may 
truly be said to have disposed of the tabulation of the Bessel functions of orders 0 and 1 
for all time. 

A description of the volume will serve to prove the truth of this statement. The Bessel 
functions Jo(x) and J:(x) are solutions of the differential equation x*y"+xy' + (x*—n*)y=0 
for n = 0 and m = 1 respectively, which are finite at x = 0. The bulk of the volume under 
review (170 out of 288 pages) is devoted to the tabulation of these functions. Ten-place 
values of Jo(x) and J:(x) and second central differences are tabulated at intervals of .001 
in the range from 0 to 16 and at intervals of .01 in the range from 16 to 25. Beyond x = 25, 
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the auxiliary functions Ao(x), Bo(x), A1(x) and Bi(x) appearing in the asymptotic expansions 
Jo(x) = Ao(x) sin x + Bo(x) cos x; Ji(x) = B,(x) sin x — A:(x) cosx are tabulated in the 
range from 25 to 1150 at intervals (ranging from .1 to 10) so chosen that second central 
differences are adequate for the maximum attainable accuracy. In the range from 1000 to 
6000 the auxiliary functions are tabulated at intervals of 100 together with ‘‘modified” 
second central differences. Beyond x = 6000, one term each in the expansions of the auxil- 
iary functions suffices to generate values to eight-place accuracy. Thus to all intents and 
purposes the volume under review disposes of the tabulation of Jo(x) and Ji(x) over the 
entire range of the argument from 0 to infinity. The situation is quite similar with the other 
Bessel functions tabulated in the volume with some exceptions to be discussed. 

The functions of the second kind Yo(x) and Yi(x) have a logarithmic singularity at the 
origin and a table of these functions would not be interpolable near the origin unless the 
interval is very small. For this reason the authors have tabulated 8-place values of the 


functions Yo(x) — = sole) Inx = Co(x) and x { Yi(x) — = J4(e) In =} = Ci(x) in the 


2 2 
range from 0 to .5 at intervals of .01. The functions Do(x) = 7; Jo(x) and Di(x) = =m 1(x), 
Tt. 


where M = log ¢, are also tabulated alongside the values of Co(x) and C(x). With the aid of 
the four auxiliary tables just described, one may compute Yo(x) = Co(x) + Do(x) log x and 
Yi(x) = Ci(x)/x + D,(x) log x. As is the practice throughout the entire volume, second 
central differences are tabulated with the values of the auxiliary functions in question. 

In addition to these functions the values of Yo(x) and Y1(x) are also given in the range 
from 0 to .5 with or without modified second central differences. Beyond x = .5 and through 
x = 25, eight-place values are given (with ordinary or modified second central differences) 
for arguments at intervals of .01. Beyond x = 25, the auxiliary functions Ao(x), Bo(x), 
A(x) and B,(x) above mentioned may be used for the computation of Yo(x) and Yi(x) by 
means of the relations Yo(x) = Bo(x) sin x — Ao(x) cos x and Yi(x) = — Ai(x) sin x 
— B(x) cos x. 

The Bessel functions I(x) and I;(x) are solutions of the differential equation x*y”’ + xy’ 
— (x? + n*)y = 0, for nm = Oand nm = 1 respectively, which remain finite at x = 0. Eight- 
place values of these functions and second central differences are given in the range from 
0 to 5 at intervals of .001. Beyond x = 5, the related functions e~*Jo(x) and e~J,(x) are 
tabulated at intervals of .01 from x = 5 to x = 10 and at intervals of .1 from x = 10 to 
x = 20. 

The functions Ko(x) and K,(x), the second fundamental solutions associated with 
Io(x) and I(x), have a logarithmic singularity at the origin and therefore for the sake of 
interpolability the functions Eo(x) = Ko(x) + Io(x) Inx and Ex(x) = x{Ki(x) — I(x) In x} 


1 
are tabulated together with the products Fo(x) = — 7 Io(x) and Fi(x) = 7; I(x) for x rang- 


ing from 6 to .5 at intervals of .01. In addition Ko(x) and Ki(x) are given for x ranging from 
0 to 5 at intervals of .01. Beyond x = 5, the values of e*Ko(x) and e*K:(x) are given along- 
side the values of e~*Jo(x) and e~*J,(x) in the same range and for the same interval as the 
last-named function. 

To round out the description of the volume under review, let it be mentioned that it 
contains the first 150 zeros of Jo(x) and J:(x) and the corresponding values of each of the 
functions at the roots of the other; the first 50 zeros of Yo(x) and Yi(x) together with the 
corresponding values of Yi(x) and Yo(x) at these zeros, a one-page table of e* and e*, a 
one-page table of the Everett coefficients of the second difference, and a 2} page table of the 
Besselian coefficients of the double second difference. 

In addition to a general introduction devoted to a description of the tables and of their 
preparation, each of the first eight major tables is preceded by a page giving the definition 
of the functions, recurrence relations, and other pertinent information. The typography of 
the volume is excellent and as far as accuracy is concerned the reviewer has the feeling 
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after reading the Introduction that the authors were fully justified in their remark that 
“There is every reason to believe that the tables are completely free from error.” 
ARNOLD N. Lowan 


180[L].—F EDERAL TELEPHONE AND Rapio Corp., Reference Data for Radio 
Engineers, New York, 1943 (fourth printing Aug. 1944), p. 197-200. 
13.8 X 21.6 cm. 


Here are 6 tables of Bessel functions: T. 1-2, Jo(x) and Ji(x), x = [0(.1)15.9; 4D]; 
T. 3-5, Jn(x), m = 2, 3,4, x = [0(.1)4.9; 4D]; T. 6, J.(p), p = 1(1)14, m = [0(.5)10; 4S]. 


n! 


181[L]—NYMTP, “Table of f,(x) = 
v. 23, 1944, p. 45-60. 17.5 X 25.3 cm. 


Under BAASMTC 2 (MTAC, p. 283) reference was made to elaborate ms. tables of 
J,(x) with the following ranges: 


n = 0(1)20, x = [0(.1)25; 10D]; n = 0(1)12, x = [0(.01)10; 8D]. 


Though these tables were made available to the authors of this table of f,(x), they have 
preferred to recompute and publish a new table, since although J,,(x) may be given to 8 or 
10 decimals, over a considerable range there are but one or two significant figures, whereas 
the function here adopted has all nine decimals significant over almost the whole range, 
particularly for m large, where the defect of the J, table is most noticeable. 

The first tables of these functions were calculated by A. WALTHER, S. GRADSTEIN & K. 
HESSENBERG and published in JAHNKE & EmpE, second ed., 1933, p. 250-258, with the 
notation 


A,(x) [J & E] = f.(x)[NYMTP], for mn = 0(1)8, x = [0(.02)9.98; 4-5D]. 


There appears no advantage in changing from the symbol A introduced by Emde, to the 
symbol f, which may even cause confusion as it is commonly needed for a general functional 
symbol. The symbol A will therefore be used in this review. ; 

This table is not only of great value for the actual numerical values thereby made 
available, but also is of theoretical interest. The functions selected for tabulation and the 
method of starting with the largest values of m and working down to small n, are clearly 
the means by which nature intended that Bessel Functions should be computed. 

The A,,(x) functions have also a valuable property perhaps not envisaged by the authors, 
opening up possibilities of double-entry tables interpolable in m as well as in x. If curves be 
plotted of An(x) against m for a fixed parameter x{1(1)10], it is seen that for m > x or there- 
abouts the curves (monotonic, asymptotic to 1 for n = ©) suggest this possibility. To test 
this, values of A,(5), m = 5(1)20, extracted from the tables and differenced, were inter- 
polated for m = 12.5, 13.5, 14.5 using 8th difference Everett (the 8th differences only affect 
a unit in the 9th (last) decimal). Compare my inefficient efforts to interpolate J,(1) between 
n = 1and 2to7D, requiring 6 differences, although fundamental values had been laboriously 
computed at interval .1 [MT AC, p. 99, where the interval erroneously given as .01 was later 
corrected, p. 132]. The three values at » = 12.5, 13.5, 14.5 satisfied the recurrence relation 
An-1(x) = An(x) — [x*Anyi(x)]/4n(m + 1); and repeated applications thereof yielded 
Az.s(5) = .08083 8726. 


‘ 3 
Now = (=) Jas(5) = (= 7 1) sin — cos 5, on making the obvious 


substitution of I'(m + 1) for m! when m is not integral. Taking sin 5 and cos 5 from 
BAASMTC, Mathem. Tables, v. 1, and computing the true value of A2.s, the result is 
.08083 87261, with which the value obtained by interpolation agrees so far as it goes. Hence 
it appears that the table can be interpolated for fractional m to nearly the full accuracy of 
the original table, provided we start from m > 2x or thereabouts. 


J,(x)”, J. Math. Phys. M.L.T., 
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So it seems that it would be possible, without excessive labour, to derive a double-entry 
table for m = 0(.1)20, x = [0(.1)10; 7D], that would be interpolable in both m and x without 
going further than second differences at most. Since A,(x) is an even function of x, it would 
be unnecessary to tabulate for x negative (compare J,(x), for fractional real m and negative 
real x, complex). But difficulties arise for A, when m is negative, since A, = © for nm a 
negative integer. Even conversion to J, is of little practical use here |J,,(.1) fluctuates from 
zero to —38 and back, between nm = — 1 and —2]. The only practical solution would appear 
to be interpolation (if required) in the region m positive, followed by the use of recurrence 
formulae. 

In equation (3), p. 46, the ” in the denominator of the last fraction should be deleted so 
as to read x?/4(m + 1)(n + 2). 


Engineering Laboratory 
University of Cambridge 


182[L].—G. N. Watson, A Treatise on the -Theory of Bessel Functions. 
Second edition, Cambridge: at the University Press, New York: The 
Macmillan Co., 1944, viii, 804 p. 18.7 X 23.2 cm. $15.00. 


This offset print of the first edition of 1922 seems to have been in the press for more than 
three years since Watson’s preface to the second edition is dated ‘‘March 31, 1941.” It is 
as follows (apart from acknowledgments of assistance) : 

“To incorporate in this work the discoveries of the last twenty years would necessitate 
the rewriting of at least chapters XII-XIX; my interest in Bessel functions, however, has 
waned since 1922, and I am consequently not prepared to undertake such a task to the 
detriment of my other activities. In the preparation of this new edition I have therefore 
limited myself to the correction of minor errors and misprints and to the emendation of a 
few assertions (such as those about the unproven character of Bourget’s hypothesis) which, 
though they may have been true in 1922, would have been definitely false had they been 
made in 1941.” 

Hence the numbers of pages in the first and second editions are the same. In the Bibliog- 
raphy, p. 753-788, the only change seems to be an addition of a title on p. 788. We note that 
equations (3), (4), p. 81 have been corrected. In MTAC, p. 307, we have already listed the 
tables in this work, p. 666-752. The errors which we there noted in T. I-II have now vanished 
but many others still remain; see MTAC, p. 296, and MTE 58, 60, where a beginning has 
been made in listing such errors. This reprint has filled a great need. RCA 


183[M].—H. W. LinpEmann, “Innenbackenbremsen,”’ Automobiltechnische 
Z., v. 46, Aug. 1943, p. 367-372; also English transl., ‘Internal shoe 
brakes,” Engineers’ Digest (Amer. ed.), v. 1, Sept. 1944, p. 558-563. 
19.7 X 27.2 cm. English ed., v. 9, Aug., 1944, p. 230-235. 

On p. 562 is a three-place table of eight integrals, from 0 to 0, of (a) sin ¢ cos ¢; (b) sin? ¢; 

(c) cos* (d) cos #; (e) sin ¢; (f) sin? ¢ cos t; (g) sin cos? t; (h) sin* for @ = 0(2°)180°. 


MATHEMATICAL TABLES—ERRATA 


References have been made to Errata in RMT 169 (BARLow, Havasu, 
WEIGEL, L. ZIMMERMANN), 170 (LEHMER), 172 (N.D.R. Comm.), 176 
(Davis), 181 (NYMTP); N26 (AcEToNn). 


55. J. M. BaTEs, ‘‘Zeros of a class of polynomials associated with Bateman’s 
k-function,” Iowa State College J. Sci., v. 12, 1938, p. 474. 


The first 5 zeros of J1/3[(2/3)x*/2] + J_1/:[(2/3)x*/] are here given as follows: 
2.338107, 4.137258, 5.520555, 6.786701, 7.944136. 
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These values should be: 
2.338107, 4.087949, 5.520560, 6.786708, 7.944134. 


Bates states that he calculated his zeros from those for the function J1/3(x) + J_1::(x), 
given by Watson, Bessel Functions, 1922 [and 1944], p. 751; but, as noted below, Watson's 
second zero is highly erroneous. 

J. C. P. Miter 


56. A. R. Forsytu, Messenger Math., s. 2, v. 50, 1921, p. 129-149. Compare 
MTAC, p. 215, IB.6. 


This paper gives exact values of certain sums extending over the roots of Bessel func- 
tions. 


P. 145, for = 73/(5-27-4096), 
read ur = 73/(5-27-4096) ; 

p. 146, for = — 37/(3*-128), = — 379/(3-5*-2048), 
read + = — 37/(3-128), 


+ = — 379/(3-5-2048). 
r=1 
To the list on p. 139 we may add the following: 


though this is not obtainable by the author’s method. 
D. H. L. 


57. K. Karas, “Tabellen fiir Besselsche Funktionen erster und zweiter 
Art mit den Parametern » = + 2/3, +1/4, +3/4,” Z. f. angew. 
Math. u. Mech., v. 16, 1936, p. 249-251. It has been already noted 
(MTAC, p. 295) that NYMTP discovered that in half of these tables, 
those for J,(x), there are 108 errors, namely: 90 of a unit in the last 
decimal place, and 18 of two or more units. These errors are as follows: 


» = 2/3 
for read x for 
1 .15011 -15012 2.1 -52905 
.23723 .23722 4.2 —.28517 
8 .54534 .54535 4.8 — .36057 
1.2 62892 .62893 5.0 — .35712 
1.6 62955 .62956 79 .27822 
1.7 61777 .61778 
y= — 2/3 
= read x for 
1 2.72981 2.72976 3.0 — 45760 
2 1.68082 1.68084 8.5 —.21775 
4 96247 -96246 8.6 — .23186 
5 76835 -76834 8.9 — .25916 
8 37691 .37690 9.3 — .25895 
2.6 — 48311 — 48310 9.6 — .23233 


read 
-52906 
— .28518 
— .36058 
— .35713 
.27823 
read 
— 45759 
—.21774 
— .23185 
—.25915 : 
— .25894 
— .23232 


| 
& 


Ue © 0-7 
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y=1/4 
for read x for read 
7143 -7144 6.9 .2548 .2547 
5803 -5804 7.0 .2681 -2680 
5383 .5384 7.2 2861 -2860 
2443 .2444 2907 -2906 
— .2189 — .2188 74 2925 -2923 
— .3836 — .3835 7.5 2913 -2910 
— .3789 — .3788 7.6 2873 .2869 
— .3436 — .3437 7.7 2805 .2800 
— .3046 — .3045 7.8 2711 .2704 
— .1345 — .1346 79 2591 -2582 
— .0684 — .0685 8.0 2449 2436 
y= —1/4 
for read = for read 
1.4310 1.4319 ZA 2742 .2743 
1.1560 1.1559 pe 2415 -2416 
— .3668 — 3669 74 2219 .2220 
— — 4226 7.5 2003 -2004 
— 4243 — 4242 7.6 1770 1771 
-3018 .3019 7.7 1522 .1523 
.3069 .3070 7.8 1263 -1264 
.2967 -2968 7.9 0994 .0995 
2868 -2869 8.0 0719 .0720 
vy = 3/4 
for read = for read 
4986 4987 —.2723 —.2722 
5988 .5989 6.1 — .1693 — .1694 
6115 .6116 6.7 0145 .0146 
6224 .6225 6.8 0447 .0449 
5888 5889 6.9 0741 .0743 
1719 1718 7.0 1024 -1025 
1274 7.3 1778 1779 
0830 .0829 7.4 1990 1991 
— .0855 — .0856 75 2180 2182 
— .1934 —.1935 7.6 2347 2348 
— .2243 — .2244 7.7 2488 2489 
— .2524 —.2525 7.8 2602 2604 
— 3439 — 3440 7.9 2689 2692 
— 3482 — .3481 8.0 2749 2752 
v= — 3/4 
Sor read x for read 
3277 
2194 .2193 3.9 —.1509 —.1508 
— 4211 — 4212 4.1 —.0719 — .0718 
— 4672 — 4673 5.3 .2941 .2942 
— .5008 — .5009 5.7 -3306 -3305 
— 4204 — .4205 6.3 -2889 .2890 
— 3943 — .3945 6.5 .2526 .2527 
— .3656 — .3657 .0998 d 
—.3007 — .3006 7.9 —.1257 — .1258 


See RMT 72, MTAC, p. 8. 


E. C. J. v. Lommet, Bayer. Akad. d. Wissen., math. natw. Abt., Abh., 


v. 15, 1886, p. 648; also G. N. Watson, A treatise on the Theory of Bessel 


Functions, 1922 and 1944, p. 744; see RMT 182 and MTAC, p. 296. 


Tables of S(2x/m)! and C(2x/m)}: 
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S(2x/x)4 C(2x/x)4 

= for read = for read 

0.2 .023788 (Watson .023720 correct) 11 90 -380392 

1.5 415348 415483 11.5 395149 .395152 

6.5 -347099 -347100 15.0 569335 -569336 

8.0 -512010 .512009 15.5 524009 -524010 
11.0 -504784 .504786 17.5 406589 406590 
11.5 447809 447810 20.5 587849 -587848 
12.0 -405810 405811 21.0 573842 .573841 
13.5 432489 432488 21.5 542266 -542265 
15.5 -598183 .598184 25.5 526896 526897 
20.5 504875 -504874 34.0 537026 .537027 
21.0 545885 .545884 34.5 504881 882 
24.0 467029 -467028 38.5 545560 545559 
28.5 573060 -573059 47.5 479313 A79311 
34.0 .557490 .557489 48.5 443930 -443929 
36.5 476871 476872 The 14 errors for x < 20 were noted 
39.5 -513690 -513689 by comparison with the tables of J. R. 
40.5 -558799 -558800 Arey, in B.A.A.S., Report, 1926, p. 
44.5 -447720 447721 274-275. 


J. W. Wrenca, Jr. 
4211 Second St., N. W. 
Washington 11, D. C. 


59. NYMTP, Tables of Circular and Hyperbolic Sines and Cosines, 1939. 
See RMT 89, p. 45f. 


The integral part in the value of sinh x corresponding to the argument 1.8190, on p. 364, 
should be 3 instead of 2. 
A. N. Lowan 


60. G. N. Watson, Bessel Functions . . . , p. 751. See RMT 182, p. 364. 


Sn, zeros of J_1/3(x) + Jiss(x) 
$2, for 5.6101956, read 5.5101956. 
J. C. P. Miter 


UNPUBLISHED MATHEMATICAL TABLES 


Reference has been made to an unpublished table in RMT 168 (TERRILL 
& SWEENY). 


29[C].—GrorGE Hastincs HEPPEL (1793-1845), Table of Logarithms, mss. 
and stereotyped plates possibly in possession of the Institute of Actuaries. 


The principal tables in these mss. consist of the logarithms, to 12D, of 100 000 to 118 600, 
stereotyped in quarto; of 131 000 to 136 636; and of 200 000 to 210 000, to 15D as far as 
202 024, and to 12D thereafter. Mr. Heppel was at various times an accountant and an 
actuary. Having decided shortly before he died that he had sufficient material to publish 
he got as far as having stereotyped plates made for the range indicated above. The accuracy 
of the tables has not been tested. The author’s son has told us that the logarithms were con- 
structed on the basis of Hutton’s logarithms of numbers to 20 places. Mainly through the 
influence of A. DEMorGan the tables and plates were purchased and placed in the hands of 
CHARLES JELLICOE, later president of the Institute of Actuaries. It would seem as if all that 
is known of the tables and their author is given in an article by Jellicoe, ““Heppel’s log- 
arithms,” Institute of Actuaries, J., v. 10, 1862, p. 82-84. 


R. C. A. 
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30[C].—Joun THomson (1782-1855), Table of Twelve-Figure Logarithms of 
Numbers from 1 to 120000, mss. presented to the Royal Astronomical 
Society in 1873 by his sister CATHERINE THOMSON. Mss. in the R.A:S. 
Library, though not listed in its printed Catalogue of 1886 or 1925. 


All that is known about these tables is contained in an ‘‘Account’”’ drawn up at the re- 
quest of the Council of the R.A.S. by J. W. L. Glaisher, and published in R.A.S., Mo. 
Notices, v. 34, 1874, p. 447-475. A few extracts from the information thus supplied are now 


The tables in about a dozen volumes include a complete 12-place table of the logarithms 
of numbers from unity to 120 005 (besides the logarithms of 744 higher numbers at irregular 
intervals, the highest being log 123 187) to 12D; first and second differences for numbers 
from 29 561 to 120 001; and logarithms of numbers from unity to 4 000 to 15D. Another 
volume containing second differences for numbers from 10 000 to 29 560 has probably been 
lost. Here are wholly original tables. 

To determine the accuracy of the mss. Glaisher read parts of them with existing tables 
giving more than 12 figures. One such table is Briggs, Arithmetica Logarithmica, London, 1924 
with 14-place logarithms of numbers from 1 to 20 000 and from 90 000 to 101 000. It was 
found that most of the ms. errors were in end-figures. By comparison with a number of other 
tables Glaisher assembled many interesting results. 

Little is known of Thomson except that he was born at Strachur, Argyleshire, and was 
the second in a family of eleven children of a farmer. He was educated at Strachur till he 
was twelve years of age, when he went to Greenock, and was placed under the care of a 
teacher of mathematics and navigation, from whom he imbibed a taste for mathematical 
studies. For many years he occupied the position of clerk to a firm in Greenock, but sub- 
sequently he commenced business as an accountant. 


31[C, D, Q].—Arithmetic, Logarithmic, Trigonometric, and Astronomical 
Tables, computed, 1848, 1869-89, by Epwarp SANG, and his daughters 
Jane Nico SANG, FLorA CHALMERS SANG, and presented in 1907 to the 
Royal Society of Edinburgh, in custody for the British Nation. 


These manuscripts, now at the R.S.E., consist of 47 volumes which were originally 
accompanied by transfer duplicates of 33 of the volumes; these latter volumes were de- 
posited in the Library of the University of Edinburgh. Of the 47 volumes, 26 (1-6, 12-19, 
28-30, 39-47) are the work of E.S., 5 (7, 32, 33, 36, 37) of J.N.S., and 16 (8-11, 20-27, 31, 
34, 35, 38) of F.C.S. The chief tables in these volumes are as follows: 

The first 10 000 Prime Numbers (up to 104 729) with their logarithms to 28D. Computed 
1848, 1865-75. 

Logarithms of Numbers 100 000-370 000, to 15D, A?. Computed 1869-1873? Also the 
last 10 figures of the logarithms of numbers from 1 00000 0000 to 1 00000 9999, and from 
1 00000 0000 to 99999 0000. Computed 1884. 

Sines for each centesimal 5 minutes of the quadrant, to 33D, A*. Computed 1877. See 
MTAC, p. 36f. 

Log Sines for every .0001 radian up to .25, and log [arc/sine], to 7D, with A. Computed 
1880-85. 

Sines and cosines for each centesimal minute in the quadrant, to 15D, A*. Computed 
1881. ; 

Log Sines for each centesimal minute in the quadrant, to 15D, 0 to 508, A; 50¢ to 100s, 
A*. Log Tangents, for each centesimal minute, to 15D, 0 to 50¢, A. Computed 1885-88; 
revised 1889. 

Sines for each centesimal minute in the quadrant, to 8D, A?. Computed 1879. 

The tables so far listed, along with their computations, are in v. 1-44. Then follow: 

Tables of Circular Segments and Mean Anomalies. Writing the equation for Kepler’s 
problem of finding the mean anomaly (m) of a planet's orbit from the eccentric anomaly (u) 
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in the form 
m=u—esinu = u — cos@sinu 


= 3{(u + 6) — sin (u + 6)} + 3{(u — sin (u — 


it appears that if u and @ are known, m may be derived from a table which gives the quantity 
4(x — sin x) with argument x. This quantity is the area of the segment of angle x* in a 
circle of unit radius. Sang first calculated a table of circular segments for each centesimal 
minute in the 400¢, to 8D, A?. The results are given in terms of the unit one grade of area, 
the sector of a circle in which the arc subtends one grade; hence the whole circle contains 
400¢ or units of area. 

Thus were derived the corresponding mean anomalies, to 8D, by a double reference, for 
each arc u = 0(1¢) 200¢, and each degree of arc of ellipticity, e = 0(1) 99«. 

Such matters are discussed, and some printed tables to 4D are given by Sang, in his 
“Nouveau calcul des mouvements elliptiques,’’ Accad. Sci., Turin, Mem., v. 32, 1880, p. 
187-199, 305-307. The table of segments is given p. 192-193; and the tables of anomalies for 
e = cos 108, cos 502, cos 90¢, p. 194-199. 

And in v. 47 “‘the Anomalies are given only to the nearest second, but the differences for 
a change of 1£ of position, and the variations for a change of 1£ in ellipticity, are filled in; 
and thus, of the three—the eccentricity, the position, the anomaly—any one may be de- 
termined frou the others.’’ The previous table to 8D is here cut to 4D, with A. 

In this same volume are also a Table of all Rational Right-angled Triangles having the 
hypotenuse less than 1000, arranged according to the magnitude of the hypotenuse; and a 
Table of Other Triangles whose sides (< 1000), and areas, are integers. 

Greater details concerning these manuscripts, and ideas developed in them, may be 
found in 

1. “Dr. Edward Sang’s logarithmic, trigonometrical, and astronomical tables,”’ R. So. 
Edinburgh, Proc., v. 28, 1908, p. 183-196. Reprinted, with a portrait photo of Sang, in 
E. M. HorssurGu, Modern Instruments and Methods of Calculation, London, 1914, p. 38-47. 

2. C. G. Knorr, “Edward Sang and his logarithmic calculations,” Napier Tercentenary 
Memorial Volume, London, 1915, p. 261-268; see also R. A. Sampson, p. 236-237. 

3. E. Sang, “On the need for decimal subdivisions in astronomy and navigation, and 
on tables requisite therefor,”” R. So. Edinburgh, Proc., v. 12, 1884, p. 533-544. 

4. E. Sang, “‘On the construction of the canon of logarithmic sines,” idem, p. 601-619; 
see also v. 9, 1878, p. 343-352. 

5. E. Sang, “Notice of fundamental tables in trigonometry and astronomy, arranged 
according to the decimal division of the quadrant,” R. So. Edinburgh, Proc., v. 16, 1889, 
p. 249-256. 

For some facts we are also indebted to Henderson (see MTAC, p. 2). 

There is a biography of Edward Sang (1805-1890) by D. B. PeEsies in R. So. Edin- 
burgh, Proc., v. 21, 1897, p. xvii—xxxii; this includes a list of his writings. Born in Kirkcaldy, 
Scotland, he was a pupil of Edward Irving (Scottish church divine and friend of Thomas 
Carlyle) until 1818, when he went to Edinburgh and entered the University. With the 
exception of two years (1841-1843) spent in Manchester New College as professor of mechan- 
ical science, and of twelve years (1843-1854) spent in Turkey, where he assisted in establish- 
ing schools of civil engineering, and in laying out railroads, he lived and worked in Edin- 
burgh as a private teacher of mathematics. For many years he was secretary of the Royal 
Scottish Society of Arts. In 1883 the University of Edinburgh conferred on him the degree 
of LL.D., and in 1884 he was elected an honorary member of the Franklin Institute in 
Philadelphia. His first published work was a large folio volume, Life Assurance and Annuity 
Tables with a copious collection of rules and examples, Edinburgh, 1841, second v. 1859. For 
interesting comments on this work, see A. DE Morcan, English Cyclopedia, ed. C. KNIGHT, 
Arts and Science section, v. 7, London, 1861, col. 1013. 

He was the first one to publish a logarithmic table beyond 108 000, A New Table of 
Seven-Place Logarithms of all numbers from 20 000 to 200 000, . . . prepared under the auspices 
of the managers of the Associated Life Insurance Offices in Scotland, London, 1871, second 
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edition, improved, Edinburgh, 1878. A set of stereotype plates for this work was given to 
the R.S.E. with Sang’s manuscripts. See J. HENDERSON, Bibliotheca Tabularum Mathe- 
maticarum, 1926, p. 135f; E. Sane, “Account of the new table of logarithms to 200 000,” 
R. So. Edinburgh, Trans., v. 26, 1872, p. 521-528. Also Specimen Pages of a table of the 
Logarithm of all Numbers up to One Million . . . shortened to nine figures from original calcu- 
lations to fifteen places of Decimals, [1874?], 26 p. 4to, not in the List of his Writings. There is 
a copy in the Brown University Library. The first page is a copy of a memorial of the 
president and council of the R.S.E. dated 19 June 1874, to the Chancellor of Her Majesty’s 
Exchequer, and the last 16 pages are devoted to communications from academies and 
scholars, dated 1872-73, commending the idea of publishing a 9-place table of logarithms of 
numbers to 10°, in three volumes; but the publication was never achieved. 

Sang “received a grant of £100 per annum frorn [the] Government as a recognition of his 
valuable scientific work; and the associated Scottish life assurance offices, feeling that some 
substantial recompense was due to him for his logarithms and actuarial tables, at a meeting 
in 1878 resolved to recommend to the offices the payment of an annuity of £100 . . . for 
the remainder of his life, which was agreed to, and subscribed by the offices.” 

R. C. A. 


MECHANICAL AIDS TO COMPUTATION 


12[Z].—L. J. Comrie, “Mechanical computing,” Appendix I, p. 462-473 
of Davip CLarkK, Plane and Geodetic Surveying, v. 2, third ed. revised by 
JaMEs GLENDENNING, London, Constable, 1943. 14 X 21.7 cm. Also as 
a pamphlet reprint. 

Some computations possible on several calculating machines are described and there 
are illustrations of four of them, namely: (1) Brunsviga, (2) 10-key Facit, (3) Marchant 
Electric, (4) Twin Marchant. Two pages are devoted to a listing and discussion of published 
tables, and some other literature. 


13[Z].—L. J. Comrtr, “Recent progress in scientific computing,’ J. Sci. 
Instruments, v. 21, Aug. 1944, p. 129-135, illustrated. 


The substance of this paper was delivered 3 July 1943 at a joint meeting of the London 
and Home Counties’ Branch of the Institute of Physics, and the London Mathematical 
Society. The headings in the paper are as follows: “Differential analysers,” ‘The training 
of computers,” ‘Numerical integration by hand,” “Finite differences,” ‘‘Direct and inverse 
interpolation,” ‘‘Double-entry interpolation,’ ‘‘Punched-card machines,” “Solution of 
simultaneous equations,” ‘“‘Mathematical tables,” “‘Short bibliography.” 


NOTES 


26. AGETON’s METHOD.—Navigators and other persons who use Hydro- 
graphic Office, Publication, no. 211 (see MTAC, p. 80-81) regularly, will 
be interested in two notes recently published on Ageton’s method of celestial 
navigation. The first is “‘The accuracy of Ageton’s method in celestial 
navigation” by SAMUEL HERRICK, Astron. So. Pacific, Publications, v. 56, 
1944, p. 149-155. Herrick points out that the warning appearing in H.O. 
211 against the use of sights for which the value of K is found to lie between 
the limits of 87° 30’ and 92° 30’ is not adequate. He shows that, when the 
tables are used in the conventional fashion and K lies outside the forbidden 
range, the maximum error in the computed altitude, 4, is six minutes of arc. 
He shows also that an error of a minute of arc in kh may occur even when K 
is quite small. He points out that the “‘caution’”’ mentioned above is unfor- 
tunately omitted from Ageton’s Manual of Celestial Navigation, New York, 
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Van Nostrand, 1942, T. II, and from three texts which reproduce the tables. 
Herrick concludes that if a maximum error of two minutes of arc in h is 
acceptable, Ageton’s method can be used without interpolation in the tables 
so long as K lies outside the range, 82° to 98°. 

The second note is ‘Increased accuracy with Ageton’s method” by 
C. H. SmiLey, Pop. Astr., v. 52, 1944, p. 379-383. Smiley points out that 
Ageton’s method could be used in many of the problems of practical as- 
tronomy if its accuracy were increased slightly. He notes that seven-place 
logarithms of sines and cosines can be used in place of Ageton’s A’s and B’s, 
otherwise keeping the computing form and rules as given by Ageton. He 
points out that the indeterminacy in the computed altitude, 4, resulting 
from the value of K being near 90° can often be avoided by interchanging 
the values of the declination and the latitude. There is a brief discussion of 
the probability that K will lie near 90°. 

C. H. SMILey 
Brown University 


27. EARLIEST TARLES OF S AND T.—In N22 A. Fletcher refers to early 
tables of the auxiliary functions S and T for arguments 0°.01(0°.01)2°.99 
given by John Newton in 1658 in his Trigonometria Britanica, and quotes a 
suggestion that these may be the first published tables. Tables embodying 
an equivalent idea, although only for arguments 0(0°.01)0°.5, were, how- 
ever, given in the previous year by William Oughtred in the second part of 
his Trigonometria,: 1657, which has a separate title-page Canones Sinuum, 
Tangentium, Secantium: et Logarithmorum pro Sinubus et Tangentibus, Lon- 
don, 1657. On p. 235 Oughtred gives quantities which he labels “SS” and 
“A” (being respectively negative and positive); these quantities, given in 
units of the seventh decimal, are respectively S — S, and T — To, in which 
Sp and 7» are the limiting values of S and T for zero argument. The end- 
figure is not completely reliable (although never more than a unit in error), 
so that it is not easy to decide whether (i) S and Sp, T and T> were indi- 
vidually rounded-off before subtraction in pairs, or whether (ii) the sub- 
tractions were first performed with more decimals and the results rounded- 
off. The former seems more likely since, if we assume this to have been done, 
Oughtred has 16 errors of +1 (i.e. errors in which his “‘S” is numerically 
too small and his ‘‘T”’ too large) and 7 of —1; with process (ii), he would 
have 5 errors of +1 and 23 of —1. 

Process (i) seems to imply that Oughtred regarded S and T as the pri- 
mary functions and performed the subtraction simply as a device for giving 
small numbers to tabulate and interpolate; it is the correct process if the 
quantities are always to be used in conjunction with the same Sp and 7» 
(i.e. for a specified unit) since the resulting value of S or T is then subject toa 
slightly smaller total rounding-off error. It seems improbable that Oughtred 
envisaged any other unit than that of 0°.01 which he used, but it is worth 
while to note that his “‘S’’ and “‘A,” if rounded-off by process (ii) above, 
would be independent of the unit used, except insofar as this determines the 
tabular arguments. 

J. C. P. MILLer 
1 The writer has seen three copies of this work, not all in agreement in all details. One 


of these belonged to Isaac Newton, whose library was recently purchased by the Pilgrim 
Trust and presented to the Library of Trinity College, Cambridge. 


| to 
the- 
0,” 
the 
lcu- 
e is 
the 
ty’s 
and 
s of 
his 
yme 
‘ing 
for 
L73 
by 
as 
ere 
ant 
hed 
don 
ical 
‘ing 
arse 
of : 
ro- 
vill 
ial 
‘ial 
56, 
0. 3 
the 
len 
rc. 
or- : 
rk, 


372 QUERIES 


28. Henry Goopwyn. A Tabular Series . . . (see MTAC, p. 22).—It 
is perhaps not generally realized that only a part of this table was ever 
published (one fifth). The complete title is: A/Tabular Series /of/Decimal 
Quotients/for all the/Proper Vulgar Fractions/of which/when in their lowest 
terms/ Neither the Numerator nor the Demoninator /is greater than /1000/Lon- 
don/Printed by W. Marchant, Ingram Court, Fenchurch Street /Published 
by J. M. Richardson, 23 Cornhill/opposite the Royal Exchange/1823. 
15 X 23.5 cm. 5 p. (introduction) + 107 p. (tables) + 5 unpaged leaves. 

In his preface Goodwyn writes that he is now publishing Part I of the 
complete tabular series, consisting of all decimals beginning with .0. Part II 
is to contain all beginning with .1, and so on to part V containing all decimals 
beginning with .4. He explains that it is unnecessary to print the remaining 
decimals from .5 to 1. since they can readily be found as the arithmetical com- 
plements of those already tabulated (because in efiect »/d=1—(d—n) /d). 
Thus ... “the number of entire quotients will amount altogether to 
152096, and there being an equal number of complementary (or easily 
calculated but not printed ones) the total will exhibit the vulgar fractions 
corresponding to 304192 decimal quotients.” He goes on to say that he is 
now printing Part I, containing all decimals beginning with 0.0, of the 
Tabular Series derived from the values given in his Table of Circles, . . . 
“but whether that table [i.e. of circles] shall be printed, or the present work 
completed, wil! depend altogether on the reception which this result of the 
Computer’s labour may meet with from the public.” 

Then follow the tabies from 1/1000 = .001 to 99/991 = .09989909, after 
which is printed “END OF PART I.” Then follow 5 unpaged leaves with 
2", m = 0(1)19; 5", 2 = 0(1)8; 275", m = 0(1)6, m = 0(1)6. 

It appears therefore that Part I of the Tabular Series was first printed, 
and then (also in 1823) the Table of Circles. It is to be presumed that, like 
most computers of tables, he found that sales did not cover the cost of 
printing, and no more of the Tabular Series was ever printed. 

In Glaisher’s articles on Goodwyn (see MTAC, p. 23) the incomplete 
state is noted. I imagine that Glaisher had a copy, but it is not in the 
“Glaisher Collection” in the Cambridge University Library, possibly be- 
cause the Library already possesses two copies. As Goodwyn’s daughter 
presented her father’s ‘complete works’ in 1831 (see MTAC, p. 100) and no 
subsequent parts accompanied the gift, it may safely be assumed that they 
do not exist. 

C. R. CosENs 
QUERIES 


11. CuBE Roots.—We have recently needed a table of cube roots of 
numbers NV = .1(.0001).2, to 6 or 7D. Is such a table available in print? 


R. F. Boyer 
Dow Chemical Co., 
Midland, Michigan 


12. ITEMs IN INTERPOLATION History.—Who originated (a) The use of 
the ‘throw-back’ to take account of higher differences neglected in interpola- 
tion, and (b) The use of interpolation formulae involving only even central 
differences, not advancing differences? Both these methods are generally 
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supposed to have come into use much more recently than 1826, but they are 
illustrated by Legendre in volume 2 of his Traité des Fonctions Elliptiques, 
Paris, 1826, paragraphs 672, 673, 676, as the following extracts show: 

672, p. 28. “Je remarquerai que lorsque les différences quatriémes 5*A 
sont assez grandes pour que les différences suivantes 5°A aient quelqu’in- 
fluence dans les interpolations, il conviendra de prendre 6*A — 765A au 
lieu de 6*A.” He states that, x being the fraction of tabular interval required 
in the interpolate (he is here using forward differences) the relevant terms 


x(x — 1)(x — 2)(x — 3) +2 ): then—‘‘comme 6°A est 


1.2.3.4 
censé trés petit par rapport a 6A, si l’on donne a x une valeur moyenne 3, le 
terme 6°A se reduira A — 


673, p. 29. Similarly a sixth difference can be “thrown back,” as we should 
now say, to the preceding fifth, by taking 565A —26*A in place of 565A, —? 

6 

676, p. 31. “Pour avoir le milieu entre deux termes consécutifs A, AI 
d'une suite dont les différences deviennent progressivement plus petites 
qu’une quantité donnée, il est bon d’avoir recours aux termes qui précédent 
et qui suivent les deux termes proposés.”’ 

Legendre is interpolating into the middle of an interval, and in modern 
notation the formula he gives is the following: 


fotfi 1.3 Act + Ai* 1.3.5 Ac® + 


being, as before, the mean value of obtained for x = }. 


24 32 246 128 
precisely the Everett formula for interpolation at 0.5 tabular interval. 
C. R. Cosens 
QUERIES—REPLIES 


12. Loc Loc TaBLeEs (Q 4, p. 131; QR 9, p. 336).—In answer to this 
question it would seem to be in order to draw attention to Slide Rules with 
log log scales, which permit ready solution of complicated problems in rais- 
ing of powers, and extracting a root besides obtaining values of natural 
logarithms, and hyperbolic functions. For example, to evaluate 


(.371)'"5-8; 2.877 = 73.7; x-378 = e* = .748; In31, 


and various hyperbolic functions of 3. 

The Keuffel & Esser Co., Hoboken, N. J., has manufactured (1) The 
log log duplex trig siide rule with trigonometric scales to represent degrees 
and minutes; (2) The log log duplex decitrig slide rule, with trigonometric 
scales to represent degrees and decimals of a degree; (3) The log log duplex 
vector rule. In Jan. 1943 these rules, 10 in. size, were respectively listed at 
$12.50, $12.50, $13.50; and the 20 in. size at $27.00, $27.00, and $31.50. 
Reproductions of the scales of these rules are given in K & E, Slide Rules and 
Calculating Instruments, New York, 1941, p. 313f-313i; and also in figs. 5—7 
of the last of the following works: 

C. N. Pickwortu, The Slide Rule: A Practical Manual, fourteenth ed., 
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Manchester, London, etc., 1916; “Slide rules with log.-log. scales,”’ p. 84-92, 
“log-log duplex slide rule,” p. 115. 

J. E. Toompson, A Manual of the Slide Rule, its History, Principle and 
Operation, New York, Van Nostrand, 1930; fifth printing, 1942 ; ‘The log-log 
scale and the log-log duplex slide rule,” p. 117-133; “Settings and typical 
problems involving the log-log scales,” p. 185-203. Thompson’s work is also 
in T. O. SLOANE, Speed and Fun with Figures, New York, Van Nostrand, 1939. 

H. O. Cooper, Slide Rule Calculations, Oxford Univ. Press, 1931; ‘‘Log- 
log scales,” p. 81-91, 118-125. 

L. M. KELts, W. F. Kern & J. R. BLAND, The Log Log Duplex Decitrig 
Slide Rule no. 4081, A Manual, New York, Keuffel & Esser, c. 1943, 103 p. 

R. W. FrENcH, Engineers’ Slide Rule, St. Louis, Mo., J. S. Swift, 1941; 
“Log log scales,” p. 66-83. 

J. J. Cuark, The Slide Rule and Logarithmic Tables including a ten-place 
table of logarithms . . ., third ed. including a table of natural trigonometric 
functions, Chicago, F. J. Drake, 1943 [c. 1941], ‘“The log-log scale,” p. 90-103. 

E. J. Hitts, A Course in the Slide Rule and Logarithms, Boston, Ginn, 
1943; “‘The log log scales,”’ p. 84-87. 

The first log log slide rules (linear and circular) were invented in 1814 
by Peter MarK RoGceEt (1779-1869) a physician in Manchester and London, 
and his paper on ‘“‘Description of a new instrument for performing mechani- 
cally the involution and evolution of numbers,’”’ R. So. London, Trans., 
v. 105, 1815, p. 9-28 + plates II-IV, led to his election, in 1815, as a fellow 
of the Royal Society. He was also a discoverer of the path for the knight 
moving on a chess-board so as to cover every square once and only once 
(Phil. Mag., v. 16, 1840, p. 306f, 498f). For a biography of Roget see R. So. 
London, Proc., v. 18, 1870, p. xxviii—xl. He was secretary of the Society for 
21 years. 

After Roget there were many others connected with various log log 
slide rules. For example J. A. Burpon (Acad. d. Sci., Paris, Comptes Rendus, 
v. 58, 1864, p. 573-576), F. BLanc (W. v. Dyck, Deutsche Math.-Ver., 
Katalog mathem. u. mathem.-phys. Modelle, Apparate u. Instrumente, Munich, 
1892, p. 145f) from which, in particular, values of sinh x and cosh x could be 
immediately read off, and W. ScHwETH, Ver. deutscher Ingenieure, Z., v. 45, 
1901, p. 567-8, 720. In his article on “Calculating machines” in v. 4, 1929, 
p. 553, of the Encyclopedia Britannica, fourteenth ed., the late D. BAXAn- 
DALL stated that ‘The log-log scale was reinvented and applied to the slide 
rule by Captain J. H. THomson in 1881, and by Prof. Joun Perry in 1902.” 
That this statement is incorrect so far as Perry is concerned, at least, may 
be noted by consulting Nature, v. 67, 1902, p. 141, which tells us that Perry 
was merely making Roget’s ideas more generally known. Nevertheless the 
“Peter and Perry” log log rule was manufactured by Albert Nestler of 
Lahr, Baden, Germany, and ‘‘Perry’s new slide rule’ by A. G. Thornton of 
Manchester, England (see Engineering Record, v. 58, Dec. 12, 1908, current 
news suppl., p. 37, and Pickworth’s Slide Rule, p. 91f). Pickworth describes 
also log log rules introduced by John Davis & Son, Ltd., of Derby, England 
(H. C. Dun op and C. S. Jackson, Slide Rule Notes, London, 1901); and 
A. W. Faber, of Stein (near Niirnberg), Germany. Thomson’s rule was ex- 
hibited in the International Exhibition of Inventions, London, 1885, with 
full recognition of its earlier discovery by Roget (Van Nostrand’s Engin. 
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Mag., v. 33, 1885, p. 517). More information about log log slide rules is 
contained in Baxandall’s Catalogue of the Collections in the Science Museum, 
South Kensington . .. Mathematics I. Calculating Machines and Instru- 
ments, London, His Majesty’s Stationery Office, 1926, p. 56-58. 

R. C. A. 


13. TABLES oF N*?? (Q5, p. 131; OR 8, p. 204, 11, p. 336).—A further 
contribution to the bibliography of these tables is American Civil Engineers’ 
Handbook, ed. by T. MERRIMAN and T. H. Wicern, fifth ed., New York, 
Wiley, 1930, p. 1312-1314. This table is for N = [.01(.01).1(.1)10(.5)20; 4S]. 

In MTAC, p. 336 a reference was given to T. 38 in King’s Handbook of 
Hydraulics, third ed. 1939, p. 103-112. In this same work, p. 117-121, is 
T. 40 “Discharge in cubic feet per second per foot of length over sharp- 
crested wiers, without velocity of approach correction by the Francis 
formula Q = 3.33H*/,” for H = [.0(.001)1.5(.01)6.99; 5S], the values Q of 
T. 40 being 3.33 times corresponding values in T. 38. A slide rule giving the 
values of the Francis’ wier formula, Q = 3.33H*/ (B — .2H) is illustrated 
in J. N. ARNOLD, Special Slide Rules, Purdue Univ., Lafayette, Indiana, 
Engineering Bulletin, v. 17, no. 5, 1933, p. 14-15. In this formula @Q is, in 
cu. ft./sec., the water discharge, with head of water H, of a rectangular 
notch wier of breadth B ft. Furthermore, the values of N*”, to a certain ac- 
curacy, could be read off at once from log log slide rules discussed in QR 12. 

& 
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211, 1. 20, omit of the first kind. 

. 215, 1. 17, end of line, for (p = 2, 4, 6), read (p = 4, 6, 8). 

217, C2 6, add , p = 1(1)18. 

. 220, last 1., for &, read 5»?. P. 225, 1. 12, for w_n(—x), read w,(—x). 

. 231, 1. —11, for odd integer, read odd-integer. 

. 234, D8A, 1. 5, for aa, read az 

240, 1. —3, in two places, for x4, read x*; 1. —5, for 10, read 8; 1. —7, before the integral, 
add: 

252, 1. 23-24, substitute the following sentence: The terms her,x, hei,x are given in 
Watson’s Bessel Functions, p. 81 and are used in Dwight 31; yerax, yeinx were added 
by J. C. P. MiLuer in the Liverpool Index. 

P. 256, E2, 1. 1, for (X/V)4, read $x(X/V)}. 

P. 257, after entry E14, add For improved forms of @ and ¢ given in nos. 1-14 we are in- 
debted to the Liverpool Index. 

271, for 1. 12 read The asymptotic forms of the ber, bei, ker, kei, functions and their 
derivatives are given in the natural form by Dwight 1; and 3,. A modified form quoted 
by Watson is 

272, 1. 2, read u = 4n(2r — 1 + 4s). 

. 285, 1. 32, for OSAKA and, read Oxaya &. 

. 287, 1. 10, for DinntK 15, read Dinnix 14. 

. 292, 1. 12, for f read J. 

295, 1. 19, for 1/3, Prescott, read 1/3; Prescott. 

300, NYMTP 8, 1. 1-3, delete Ko(x) = Eo(x) . . . x = 0(.001).03; and. 

. 330, 1. 11, for m = 20, read n = 20. 

. 333, 1. 16, and 17, editorial slips for which the author was not responsible: for 0(0°.001)3°, 

read 0°.01(0°.01)2°.99, and for 0°.001, read 0°.01. 
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